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ABSTRACT 

We investigate the [O II] emission-line as a star formation rate (SFR) indicator using integrated 
spectra of 97 galaxies from the Nearby Field Galaxies Survey (NFGS) . The sample includes all Hubble 
types and contains SFRs ranging from 0.01 to 100 Mq yr~^. We compare the Kennicutt [O II] and Ha 
SFR calibrations and show that there are two significant effects which produce disagreement between 
SFR([0 II]) and SFR(Ha): reddening and metallicity. Differences in the ionization state of the ISM do 
not contribute significantly to the observed difference between SFR([0 II]) and SFR(Hq;) for the NFGS 
galaxies with metallicities log(0/H) + 12> 8.5. The Kennicutt [O II]-SFR relation assumes a typical 
reddening for nearby galaxies; in practice, the reddening differs significantly from sample to sample. 
We derive a new SFR([0 II]) calibration which does not contain a reddening assumption. Our new 
SFR([0 II]) calibration also provides an optional correction for metallicity. Our SFRs derived from 
[O II] agree with those derived from Ha to within 0.03-0.05 dex. We show that the reddening, E(i? — V), 
increases with intrinsic (i.e. reddening corrected) [O II] luminosity for the NFGS sample. We apply our 
SFR([0 II]) calibration with metallicity correction to two samples: high-redshift 0.8 < z < 1.6 galaxies 
from the NICMOS Ha survey, and 0.5 < z < 1.1 galaxies from the Canada- France Redshift Survey. The 
SFR([0 II]) and SFR(Ha) for these samples agree to within the scatter observed for the NFGS sample, 
indicating that our SFR([0 II]) relation can be applied to both local and high-z galaxies. Finally, we 
apply our SFR( [O II] ) to estimates of the cosmic star formation history. After reddening and metallicity 
corrections, the star formation rate densities derived from [O II] and Ha agree to within 30%. 

Subject headings: galaxies:starburst-galaxies:abundances-galaxies:fundamental 
parameters-galaxies:high-redshift 

1. INTRODUCTION 

Observing the star formation rate since the earliest times in the Universe is crucial to our understanding of the formation 
and evolution of galaxies. The star formation rate indicators developed four decades ago provided a first quantitative 
measure of the global star formation in galaxies (Tinsley 1968, 1972; Searle, Sargent, &; Bagnuolo 1973). These indicators 
were based on stellar population synthesis models of galaxy colors. More recent and precise star formation rate indicators 
rely on optical emission-lines, UV continuum, radio, and infrared fluxes (e.g., Kennicutt & Kent 1983; Donas & Deharveng 
1984; Rieke & Lebofsky 1978). These indicators, applied to nearby samples, provide insight into the properties of galaxies 
along the Hubble sequence (see Kennicutt 1998, for a review). 

The advent of large spectroscopic surveys enabled significant progress in our understanding of global galaxy evolution as 
a function of redshift. Lilly et al. (1995) studied the cosmic evolution of the field galaxy population to a redshift of z ~ 1 
using the Canada- France Redshift Survey (CFRS). They showed that the field galaxy population evolves and that this 
evolution is strongly related to galaxy color. Ellis et al. (1996) confirmed this observation using Autofib redshift survey 
data over a similar redshift range. Ellis et al. concluded that the steepening of the luminosity function with look-back 
time is a direct consequence of the increasing space density of blue star forming galaxies at moderate redshift. 

Deep surveys like the Hubble Deep Fields allowed the study of star formation history of galaxies over an even wider 
redshift range. Madau et al. (1996) estimated the star formation history of galaxies between 2 = and z = 4. Using the 
Hubble Deep Fields and UV surveys from Lilly et al. (1996), Madau et al. argued that the peak star formation rate occurs 
at redshifts from z = 1.3 — 2.7. Many large, deep spectroscopic surveys carried out recently have sparked an explosion of 
research into the star formation history of the universe (for example. Hammer et al. 1997; Rowan- Robinson 2001; Cole et 
al. 2001; Baldry et al. 2002; Lanzetta et al. 2002; Rosa- Gonzalez, Terlevich, & Terlevich 2002; Tresse et al. 2002; Hippelein 
et al. 2003). 

Cosmic star formation history studies over a large redshift range require the use of different star formation rate indicators. 
Unfortunately, there are significant discrepancies among SFR estimates made using different indicators. To obtain a more 
reliable view of the cosmic star formation history, it is essential to gain a detailed understanding of and to reach agreement 
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among the star formation indicators at multiple wavelengths . 

The hydrogen Balmer line Ha 6563A is currently the most reliable tracer of star formation, provided Ha can be 
corrected for reddening. In the ionization-bounded nebulae of H II regions and star-forming galaxies, the Balmer emission 
line luminosity scales directly with the total ionizing flux of the embedded stars. For many years there was an apparent 
disagreement between the star formation rate derived from Ha and those derived at other wavelengths, including the 
far- infrared (FIR) . Correction of Ha for stellar absorption and reddening brings the Ha and FIR SFRs into agreement for 
active star-forming galaxies (e.g., Rosa- Gonzalez, Terlevich, & Terlevich 2002; Chariot et al. 2002; Dopita et al. 2002) and 
for the normal star- forming galaxies of all Hubble types in the Nearby Field Galaxy Survey (Kewley et al. 2002; Jansen 
et al. 2000a,b). 

Although Ha provides a useful SFR indicator for nearby galaxies, it is not easily observable for more distant galaxies. Ha 
redshifts out of the visible band for z > 0.4. An alternative diagnostic for the z ~ 0.4 — 1.5 range is the [O II] AA3726, 3729 
doublet. Several authors have calibrated the [O II] star formation rate (e.g., Gallagher et al. 1989; Kennicutt 1998; 
Rosa-Gonzalez, Terlevich, & Terlevich 2002). Unfortunately, the [O II] emission-line is plagued by problems including 
reddening and abundance dependence (Jansen, Franx, & Fabricant 2001; Chariot et al. 2002). Most previous comparisons 
of [O II] with other SFR indicators were based on spectra which lack sufficient spatial coverage, signal-to-noise, and/or 
wavelength coverage to make a detailed correction for reddening (eg., Kennicutt 1983; Hopkins et al. 2001; Chariot et al. 
2002; Buat et al. 2002). For example. Chariot et al. (2002) had to assume an 'average' attenuation Av = 1 for the galaxies 
in the Stromlo-APM survey. They found significant discrepancies between the Ha and [O II] SFRs which they attributed 
to variations in the effective gas parameters (ionization, metallicity, and dust content) of the galaxies. Similarly, Teplitz 
et al. (2003) showed that there is a large disagreement among the SFR estimates based on Ha and those based on the 
[O II] emission line. 

At temperatures typical of star-forming regions (10000-20000 K), the excitation energy between the two upper D levels 
for [O II] and the lower S level is roughly the thermal electron energy kT . The [O II] doublet is therefore closely linked 
to the electron temperature and consequently abundance. In fact, the [O II] emission-line doublet enters most reliable 
optical abundance diagnostics developed over the last two decades (eg., Pagel, Edmunds & Smith 1980; Edmunds & Pagel 
1984; Dopita & Evans 1986; Torres-Peimbert, Peimbert & Fierro 1989; Skillman, Kennicutt & Hodge 1989; McGaugh 
1991; Zaritsky, Kennicutt & Huchra 1994; Chariot & Longhetti 2001; Kewley & Dopita 2002). Most of these abundance 
diagnostics use the intensity ratio ([O II] A3727 + [O III] AA4959,5007)/H/3, commonly known as R23. Jansen, Franx, & 
Fabricant (2001) showed that the [O II]/Ha ratio is strongly correlated with the R23 ratio for the Nearby Field Galaxies 
Survey (NFGS). Jansen et al. concluded that [O II] is affected by metallicity and that H/3 is a significantly better tracer of 
star formation when detected at a sufiicient signal-to- noise (S/N) and spectral resolution to correct for underlying stellar 
absorption. None of the SFR calibrations so far take abundance into account, largely because a set of high S/N integrated 
(global) spectra for galaxies spanning a large range of SFRs is required to derive a reliable calibration. Such a sample has 
not been available. 

Here, we investigate [O II] as a star formation rate diagnostic using integrated (global) spectra for the NFGS. Our 

spectra have the advantage that: (1) they contain both [O II] and Ha, (2) we can correct Ha and H/3 for underlying 
stellar absorption, (3) we can measure the reddening using the Balmer Decrement, and (4) we can resolve Ha and [N II]. 
Furthermore, we can calculate global galaxy abundances from the common diagnostic R23 to investigate the relationship 
between SFR([0 II]) and abundance. 

We describe the sample selection and optical spectra in Section 2. The commonly-used SFR indicators are discussed in 
Section 3. In Section 4, we explore the discrepancy between the SFR([0 II]) and SFR(Ha) and derive a new SFR([0 II]) 
calibration which takes abundance and reddening into account. In Section 5, we use theoretical population synthesis and 
photoionization models to investigate the theoretical dependence of SFR([0 II]) on the ionized gas properties and we 
derive a theoretical SFR([0 II]) diagnostic. In Section 6, we discuss the use of [O II] as a SFR indicator in more distant 
samples and apply our SFR([0 II]) calibrations to the 0.8 < 2: < 1.6 sample of Hicks et al. (2002) and the 0.5 < 2 < 1.1 
sample of Tresse et al. (2002). Based on these findings, in Section 7, we investigate the implications for cosmic star 
formation history studies which use [O II] as a SFR indicator. Throughout this paper, we adopt the fiat A-dominated 
cosmology as measured by the WMAP experiment (/i = 0.72, = 0.29; Spergel et al. 2003). 

2. SAMPLE SELECTION AND SPECTROPHOTOMETRY 

The NFGS is ideal for investigating star formation rates because it is an objectively selected sample for which integrated 
spectra are available. Jansen et al. (2000a) provide a detailed discussion of the NFGS sample selection. Bricfiy, Jansen 
et al. selected 198 nearby galaxies in an objective (unbiased) manner by sorting the CfA catalog into 1 mag-wide bins 
of Mz- Within each bin, the sample was sorted according to their CfAl morphological type. To avoid a strict diameter 
limit, which might introduce a bias against the inclusion of low surface brightness galax;ies in the sample, Jansen et al. 
used a radial velocity Hmit, VLG(kms^^) > iq-O-i^-O-^Mz respect to the Local Group standard of rest). To avoid a 

sampling bias favoring a cluster population, they excluded galaxies in the direction of the Virgo Cluster. Finally, Jansen 
et al. selected every Nth galaxy in each bin to approximate the local galaxy luminosity function (e.g., Marzke, Huchra, 
& Geller 1994). The final 198-galaxy sample represents the full range in Hubble type and absolute magnitude present in 
the CfAl galaxy survey (Davis & Peebles 1983; Huchra et al. 1983). 

Both integrated and nuclear spectrophotometry are available for almost all galaxies in the NFGS sample, including 
integrated Ha, H/3, and [O II] fiuxes (Jansen et al. 2000b). The integrated spectra typically cover 82±7% of each 
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galaxy. We have calibrated the integrated fluxes to absolute fluxes by careful comparison with B-band surface photometry 
(described in Kewley et al. 2002, ; hereafter Paper I). The Ha and H/3 emission-line fluxes were carefully corrected for 
underlying stellar absorption as described in Paper I. 

A total of 116 galaxies in the NFGS have spectra with measurable Ha, H/3, [N II] A6584 fluxes, [O II] A3727, and 
[O III] A5GG7 emission lines. Due to low S/N ratios in the [O III] A4959 emission-line, we used the theoretical ratio 
[O III] A5007/[O III] A4959 - 3 to calculate the [O III] A4959 flux. 

The NFGS emission-line fluxes have been corrected for Galactic extinction by two methods: (1) using the HI maps of 
Burnstein & Heiles (1984), listed in the Third Reference Catalogue of Bright Galaxies (de Vaucouleurs et al. 1991), and 
(2) using the COBE and IRAS maps (plus the Leiden-Dwingeloo maps of HI emission) of Schlegel, Finkbeiner Davis 
(1998). The average Galactic extinction is E{B - F)=G.014 ± 0.003 (method 1) or E{B - y)=0.016 ± 0.003 (method 2). 

We corrected the emission line fluxes for reddening using the Balmer decrement and the Gardelli, Clayton, &Mathis 
(1989) (CCM) reddening curve. We assumed an Ry = Av/E(B - V) = 3.1 and an intrinsic Ha/H/J ratio of 2.85 (the 
Balmer decrement for case B recombination at T= W^K and rig ~ 10^ — lO^cm"^; Osterbrock 1989). After underlying 
Balmer absorption was removed, ten galaxies have Balmer decrements less than 2.85. A Balmer decrement less than 
2.85 results from a combination of: (1) intrinsically low reddening, (2) errors in the stellar absorption correction, and (3) 
errors in the line flux calibration and measurement. Errors in the stellar absorption correction and flux calibration are 
discussed in detail in Paper I, and are ~12-17% on average, with a maximum error of ~ 30%. For the S/N of our data, the 
lowest E(B-V) measurable is 0.02. We therefore assign these ten galaxies an upper limit of E(B-V)< 0.02. The difference 
between applying a reddening correction with an E(B-V) of 0.02 and 0.00 is minimal: an E(B-V) of 0.02 corresponds to 
an attenuation factor of 1.04 at Ha and 1.09 at [O II] using the CCM curve. 

To rule out the presence of AGN in the NFGS sample, we used the theoretical optical classiflcation scheme developed 
by Kewley et al. (2001a). The optical diagnostic diagrams indicate that the global spectra of 97/116 NFGS galaxies are 
dominated by star formation. These 97 galaxies (Table 1) constitute the sample wc analyse here. The spectra of the 
remaining 19 galaxies are either dominated by AGN (5/19) or are "ambiguous" galaxies (14/19). Ambiguous galaxies 
have line ratios that indicate the presence of an AGN in one or two out of the three optical diagnostic diagrams. Because 
these galaxies are likely to contain both starburst and AGN activity (see e.g., Kewley et al. 2001a; Hill et al. 1999), we 
do not include them in the following analysis. 

3. THE K98 [oil] AND Ho; SFR INDICATORS 

The development of SFR calibrations has been an intense topic of research for more than three decades (see Kennicutt 
1998, for a review) (hereafter K98). We start with the K98 SFR relations for [O II] and Ha because these relations are 
applied in many current SFR studies. 

The K98 Ha SFR calibration is derived from evolutionary synthesis models that assume solar metallicity and no dust. 
K98 assumed that the total integrated stellar luminosity shortward of the Lyman limit is re-emitted in the nebular emission 
lines. The K98 relation between Hq luminosity and SFR is: 

SFR(MQyr-i) = 7.9 x lO"''^ L(Ha) (ergss'^) (1) 

where L(Ha) denotes the intrinsic Ha 6563A luminosity. Paper I shows that, once the SFR(Ha) is corrected for underlying 
Balmer absorption and reddening, the mean SFR derived from Ha agrees with the mean SFR derived from the far-infrared 

luminosity to within 10% (Kewley ct al. 2002). 

The [O II] SFR calibration is much less straightforward. The K98 [O II] calibration is: 

SFR(MQyr-i) = (1.4 ± 0.4) x lO-''^ L([OII]) (ergs s"^) (2) 

K98 derived this calibration from the K98 SFR(Ha) relation (equation 1) and two previous [O II] calibrations by Gallagher 

et al. (1989, ; hereafter G89) and Kennicutt (1992, ; hereafter K92). The error estimate in equation (2) reflects the 
difference between the G89 and K92 samples. There are a number of important points to consider regarding the K92 and 
G89 calibrations: 

1. The K92 calibration is based on the [O II]/(Ha-|-[N II]) ratio. K92 assumes an average value for [N II]/Ha of 0.5 
because Ha and [N II] are blended for many galaxies in the K92 sample. Recent higher resolution spectroscopy and 
theoretical photoionization models show that the average [N II]/Ha ratio is around 0.5 ± 0.2 for most optical and infrared 
selected samples with metallicities exceeding 0.5xsolar. For metallicities below 0.5xsolar, the [N II]/Ha ratio may be as 
low as 0.01 (Kewley et al. 2001b). The mean [N II]/Ha ratio for any particular sample depends on the sample selection 
criteria and on the diagnostic used to remove galaxies containing AGN from the sample. For the NFGS, the [N II]/Ha 
ratio ranges between 0.03 and 0.50 (Jansen et al. 2000b) with a mean value of 0.27 ± 0.01, significantly different from the 
K92 value of 0.5. 

2. The K92 calibration is derived by starting with an Ha SFR calibration from population synthesis models assuming 
no dust. K92 uses the average [O II] /Ha ratio for their sample to convert the Ha SFR calibration into an SFR calibration 
based on [O II]. The average [O II] /Ha ratio used by K92 is uncorrected for reddening. Such reddening corrections were 
difficult to make at the time: lower resolution and lower S/N spectra limited the ability to obtain a reliable estimate 
for the Balmer decrement. The K92 [O II] SFR calibration therefore assumes an average reddening for the sample. The 
G89 [O II] SFR calibration is derived in a similar manner, but with an uncorrected average [O II] /H/3. The effect of 
reddening is less severe for the [O II] /H/3 ratio than for the [O II] /Ha ratio. (Note that the K98 [O II] SFR indicator 
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requires correction for reddening at Ha rather than at [O II] because the Ha SFR caUbrations are calculated from stellar 
population synthesis models assuming no dust.) 

The problem with applying the K92 or G89 [O II] indicators to individual galaxies or to samples of galaxies is that the 
reddening between [O II] and Ha may not be the same as the average reddening for either of the K92 or G89 samples. 

Indeed, Aragon-Salamanca et al. (2003) showed that prior to reddening correction, there is a significant difference between 
the average [O II] /Ha ratio for the NFGS sample and for galaxies in the Ha selected Universidad Complutcnse de Madrid 
(UCM) Survey. 

In Figure (la), we show the difference between the [O II]/Ha ratio for the NFGS and the K92 high resolution sample, 
prior to reddening correction. The mean [O II] /Ha for the K92 sample is 0.57 ± 0.06 compared to 0.73 ± 0.03 for the 
NFGS. However, after correction for reddening, the [O II]/Ha ratio difference disappears (Figure lb). The mean [O II]/Ha 
for the K92 sample after reddening correction is 1.1 ± 0.1, compared to 1.2 ± 0.3 for the NFGS. Aragon-Salamanca et al. 
(2003) found a similar agreement between the mean [O II] /Ha for the NFGS and UCM samples after reddening correction. 
These results suggest that the [O II] SFR indicator can be recalibrated in a reddening-independent manner. The SFR 
[O II] indicator would then be applicable in an unbiased way to a wider range of samples. 

3. Hidden in the [O II] SFR indicator may be errors resulting from stellar absorption undcaiying the Ha and H/3 
emission-lines. Balmer absorption is difficult to measure reliably with lower S/N, lower resolution data. This problem 
does not affect the K98 Ha SFR indicator; it assumes that the Ha emission- line is corrected for underlying stellar 
absorption. This problem does however affect the [O II] /Ha ratio derived for the K92 sample and the [O II]/H/3 ratio for 
the G89 sample. Underlying stellar absorption reduces the flux of the Ha or H/3 emission-line; without correction, the 
[O II] /Ha ratio is overestimated. 

4. Differences in metallicity also hinder the calibration of [O II]. K98 states that metal abundances have a relatively small 
effect on the [O II] calibration over most of the abundance range of interest for the K92 galaxies: (0.05 Zq < Z < 1 Zq). 
However, Jansen, Franx, & Fabricant (2001) find that there is a correlation between [O Il]/Ha and the oxygen-abundance 
sensitive line ratio R23. Chariot et al. (2002) observe a similar dependence. 

In summary, the SFR([0 II]) estimated using K98 for any individual galaxy may not provide the true SFR because of 
differences in the reddening, Balmer absorption, the [N II] /Ha ratio, ionization properties and metallicity of the galaxy 
compared to the average of the K92 sample. For an entire sample, the combination of these effects could result in an 
increase in the error (scatter) in the SFR([0 II]) relation, and possibly systematic shifts, depending on the sample selection 
criteria. 

In Figure (2), we compare the SFR(Ha) with the SFR([0 II]) derived using the K98 calibrations. We corrected the 
[O II] flux for reddening at Ha as required by K98. We flt a straight line to the logarithm of the SFRs using linear 
least-squares minimization that includes error estimates for both variables. We assumed errors of ~30% as in Kewley et 
al. (2002). The resulting fit (dotted line in Figure 2) has the form: 



Figure 2 shows that the K98 SFR(Ha) caUbration predicts a lower SFR than the K98 SFR([0 II]) calibration for SFRs 
below 1 M0/yr, but a larger SFR estimate for SFRs above 1 M0/yr. The rms dispersion around the line of best-fit in 
Figure 2 is 0.11 in the log. We will investigate the difference in slope and the relatively large scatter in Section 4. 

Other SFR([0 II]) calibrations have been derived in a manner similar to K98. Hippelein et al. (2003) provided an [O II] 



log[SFR([OII])] = (0.83 ± 0.02) log[SFR(Ha)] + (0.01 ± 0.02) 
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Fig. 1. — Histograms showing the [O II] /Ha distribution for the NFGS and K92 samples (a) prior to reddening correction, and (b) after 
reddening correction. The vertical lines at the top indicate the mean [O II] /Ha for each sample. Prior to reddening correction the means are: 
0.57 ± 0.06 (K92) and 0.73 ± 0.03 (NFGS). After reddening correction, the means are: 1.1 ± 0.1 (K92) and 1.2 ± 0.03 (NFGS). 
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Fig. 2. — Bottom panel: A comparison between the Ha and [O II] SFRs based on the K98 cahbrations (equations 1 and 2) Ha has been 
corrected for reddening using the Balmer decrement. We corrected the [O II] flux for reddening at Ho as required by K98. The legend 
indicates the Hubble type. Top panel: The Ho SFR versus the ratio of the two SFRs from Figure (a) : SFR([0 II])/SFR(Ha). In both 
panels, the solid line shows where the data would lie if both SFR indicators agreed (y=x) and the dotted line shows the least squares fit to 
the data. We assume errors of ~ 30% for SFR(Ha) and SFR([0 II]), as in Kewley et al. (2002). 
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SFR based on extinction-corrected [O II] /Ha measurements but did not correct for Balmer absorption. Rosa-Gonzalez, 
Terlevich, & Terlevich (2002), however, did correct for reddening and underlying Balmer absorption. Gallagher et al. 
(1989) and Cowie et al. (1997) used the [O II]/Ha ratio to obtain the SFR([0 II]) for different samples. None of these 
methods, however, take abundance into account. 

4. DERIVATION OF THE NEW SFR([0 II]) INDICATOR 

As mentioned earlier, in the NFGS spectra, (1) the [N II] and Ha lines are cleanly separated, (2) reddening can be 
estimated from the Balmer dccrenicint, and (3) the; stellar absorption under Ha can be measured from the H/3 emission- line 
profile. Furthermore, theoretical strong-line abundance diagnostics now enable the reliable determination of abundances 
from a wide variety of available emission-lines (eg., McGaugh 1991; Zaritsky, Kennicutt & Huchra 1994; Chariot et al. 
2002; Kewley & Dopita 2002). These diagnostics allow us to derive a new SFR([0 II]) calibration which includes an 
explicit correction for abundance. 

4.1. Reddening Correction 

In Figure (3a) we show the relationship between [O II] /Ha and the reddening E{B — V) derived from the Balmer 
decrement. The [O II]/Ha ratio is uncorrected for reddening. The Spearman Rank correlation coefficient is -0.80. The 
two-sided probability of obtaining a value of -0.80 by chance is almost zero 1.1 x 10~^^), indicating a strong correlation 
between [O II] /Ha and E{B-V). 

As described in Section 2, we correct the [O II]/Ha ratio for reddening using the CCM reddening curve. Figure (3b) 
shows the relationship between [O II] /Ha and E{B — V) after reddening correction. The Spearman Rank correlation 
correlation coefficient is -0.02. The probability of obtaining a value of -0.02 by chance is 83%, indicating that reddening 
correction by the CCM method removes the dependence of [O II] /Ha on E(i3 — \^). 

The mean [O II]/Ha for the NFGS sample after reddening correction is 1.2 ±0.3. If we apply this factor to equation (1), 
we obtain: 

SFR([OII])(MQyr-i) = (6.58 ± 1.65) x lO-'*^ L([OII]) (ergss"^) (4) 

where L([0 II]) should be corrected for reddening at [O II]. Note that in contrast to the K98 SFR([0 II]) equation (2), 
our equation (4) makes no assumption about the typical reddening. 

Figure 4 shows the SFR([0 II]) derived with equation (4) compared to SFR(Ha). The line of best-fit to the data is 
now: 

log[SFR([OII])] = (0.97 ± 0.02) log[SFR(Ha)] + (-0.03 ± 0.02) (5) 

Clearly the difference in reddening between [O II] and Ha for the NFGS compared to the K92 sample is responsible for 

the departure of the slope from unity in Figure 2. The rms scatter of the data about the fit is slightly smaller, 0.08 dex. 
In Sections 4.2 and 4.3, we explore the remaining sources of scatter between SFR([0 II]) and SFR(Ha). 

4.2. [O II]/Ha and the Ionization Parameter 
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Fig. 3. — The reddening E(_B — V) derived from the Balmer decrement compared witii the [O II]/Hq: ratio; (a) prior to reddening correction, 
and (b) after reddening correction. All fluxes have been corrected for Galactic extinction. The estimated errors axe ~ 17% and ~ 23% for 
the uncorrected and corrected [O II]/Ha ratios respectively. The estimated error in E(B — V) is ±0.04. Symbols are as in Figure 2. 
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Fig. 4. — Bottom panel: A comparison between the SFR(Ha) and our SFR([0 II]). SFR([0 II]) has been derived using our new [O II]/Ha 
ratio as given in equation (4). Both [O II] and Ha have been corrected for reddening using the Balrner decrement. The [O II] luminosity 
has been corrected for reddening at [O II], as required by equation (4). The soUd Une is y=x. The dotted Une shows the least squares 
fit to the data. The legend indicates the Hubble type. Top Panel: The Ha SFR versus the ratio of the two SFRs from the bottom panel: 
SFR([0 II])/SFR(Ha). The solid line shows where the data would lie if both SFR indicators agreed. Estimated errors are ~ 30% for SFR(Ha) 
and ~ 35% for SFR([0 II]). The dotted curve corresponds to the least squares fit. 
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There is some concern (eg K98) that the SFR([0 II]) cahbration is less precise than the SFR(Hq;) cahbration because 
[O II] is sensitive to the excitation state of the gas. For example, the excitation of [O II] is particularly high in the diffuse 
gas in starburst galaxies (e.g., Martin 1997). The ionization parameter is a measure of the excitation of the gas, and is 
defined as 



5 = 



(6) 



where 5ho is the ionizing photon flux per unit area, and n is the local number density of hydrogen atoms. The ionization 
parameter q, can be physically interpreted as the maximum velocity of an ionization front driven by the local radiation 
field. Dividing by the speed of light gives the more commonly used dimcnsionless ionization parameter; U = q/c. 

If the [O II] SFR calibration depends upon the ionization parameter, then we expect to observe this dependence in the 
[O II] /Ha ratio. In Figure 5, we plot the ionization-parameter sensitive ratio [O IIl]/[0 II] versus [O II] /Ha. The Spearman 
Rank correlation coefficient is 0.11. The two-sided probability of obtaining this value by chance is 30%, indicating that 
there is no statistically significant dependence of [O II] /Ha on the ionization parameter as traced by [O III]/[0 II]. Our 
local sample covers a small range in ionization parameter (1 x 10^ - 3 x 10^ cm/s; Dopita et al. 2001). The majority of the 
oxygen emission in the NFGS is likely to result from the 0+ species: the [O I] A6300 emission is weak or immeasurable, 
and the majority of the NFGS (72%) have [O IIl]/[0 II] ratios less than 0.5, with a mean([0 III]/[0 II])~ 0.38 ± 0.03. 

Note that the NFGS is representative of galaxies in the local universe. Samples which have not been objectively 
selected, and perhaps those at high redshifts could exhibit different ionization properties from those observed in the 
NFGS. In particular, active starburst galaxies and blue compact galaxies may contain radiation fields characterized by 
larger ionization parameters than observed for the NFGS (e.g., Martin 1997). For example, the K92 local "high resolution" 
sample (excluding the galaxies known to contain AGN) has a mean [O III]/[0 II] ratio of 0.5 ± 0.2, compared to the mean 
NFGS [O III]/[0 II] ratio of 0.38 ± 0.03. The larger K92 mean [O III]/[0 II] is caused by one galaxy in the K92 sample 
that has an extremely large [O III]/[0 II] ratio of 4.57, a factor of 8 times larger than any other galaxy in the K92 sample. 
If this outlying galaxy is removed, the average [O III]/[0 II] ratio is much lower: [O IIl]/[0 Il]= 0.28 ± 0.03. Clearly 
the [O III]/[0 II] ratios can vary significantly from galaxy to galaxy. In addition, the range in ionization parameter and 
metallicity covered by a particular sample may be infiuenced by the sample selection criterion. Lilly, Carollo, & Stockton 
(2003) observed the [O in]/[0 II] ratio for 66 galaxies with redshifts 0.47 < z < 0.92. They find that the [O III]/[0 II] ratio 
(uncorrected for reddening) in these galaxies is ~ 0.1 — 1.3. Lilly et al. note that if an average reddening of E{B — V)^ 0.2 
is applied to the CFRS sample, then the range in [O IIl]/[0 II] for the CFRS sample is similar to the range observed in the 
NFGS. However, the [O III]/[0 II] ratio is much higher in the five Lyman break galaxies (z ^ 3) observed by (Pcttini et al. 
2001) with [O II] and [O III] line fluxes. These Lyman break galaxies all have 1 <[0 III]/[0 II]< 10. The dominant process 
affecting the [O II] /Ha ratio in such galaxies may be ionization parameter rather than abundance because relatively large 
amounts of oxygen may exist in [O III] AA4959,5007 and higher levels of excitation. 

4.3. [O n]/Ha and the Oxygen Abundance 

We now investigate the dependence of the [O II] /Ha ratio on the oxygen abundance. The oxygen abundance is 
ideally measured directly from the ionic abundances obtained from a determination of the electron temperature of the 
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galaxy. An appropriate correction factor accounts for the unseen stages of ionization. The electron temperature can be 
determined from the ratio of the auroral line [O III] A4363 to a lower excitation line such as [O III] A5007. In practice, 
however, [O III] A4363 is very weak in metal-poor galaxies, and is not observed in higher metallicity galaxies. In addition, 
[O III] A4363 may be subject to systematic errors when using global spectra: Kobulnicky, Kennicutt & Pizagno (1999) 
found that for low metallicity galaxies, the [O III] A4363 diagnostic systematically underestimates the global oxygen 
abundance. 

Without a reliable electron temperature diagnostic, global abundance determinations are dependent on the measurement 
of the ratios of strong emission-lines. The most commonly-used ratio is ([O II] A3727 -|- [O III] AA4959,5007) /H/3 (otherwise 
known as R23), first proposed by Pagel et al. (1979). 

The logic for the use of this ratio is that it provides an estimate of the total cooling due to oxygen. Because oxygen 
is one of the principle nebular coolants, the R23 ratio should be sensitive to the oxygen abundance. One of the caveats, 
however, with using R23 is that it is double valued: at low abundance, the intensity of the oxygen lines scales roughly 
with the chemical abundance; at high abundance the nebular cooling becomes dominated by the infrared fine structure 
lines and the electron temperature (and therefore R23) decreases. Detailed theoretical model fits to H II regions have 
been used to develop a number of calibrations of R23 with abundance (see e.g., Kewley & Dopita 2002, for a review). 
Calibrations of R23 produce oxygen abundances which are generally comparable in accuracy to direct methods relying on 
the measurement of nebular temperature, at least in the cases where these direct methods are available for comparison 
(McGaugh 1991). 

Because different abundance diagnostics can have systematic problems, we applied four independent abundance diag- 
nostics; (1) the Kewley & Dopita (2002) [N II]/[0 II] diagnostic (hereafter KD02), (2) the McGaugh (1991) R23 diagnostic 
(hereafter M91), (3) the Zaritsky, Kennicutt & Huchra (1994) R23 diagnostic (hereafter Z94), and (4) the Chariot & 
Longhetti (2001) "case F" diagnostic (hereafter COl). 

The KD02 [N II] /[O II] calibration is based on a combination of stellar population synthesis and detailed photoionization 
models. The [N II]/[0 II] ratio is sensitive to abundance for log(0/H)-|-12> 8.5 for two reasons: (1) [N II] is predominantly 
a secondary element for log(0/H) -|- 12> 8.5, and therefore [N II] is a stronger function of metallicity than [O II], (2) at 
high metallicity, the lower electron temperature decreases the number of coUisional excitations of the [O II] lines. For 
log(0/H) -I- 12< 8.5, [N II]/[0 II] is less sensitive to abundance, and is only useful for providing an initial guess to a more 
sensitive abundance diagnostic. There are 17/97 galaxies with log(0/H) + 12< 8.5 (log([N II]/[0 II])< —1.1). Four galaxies 
have very low [N II]/[0 II] ratios (log([N II]/[0 II])< -1.43). Only galaxies with low abundances (7.5 <log(0/H)-M2< 8.2) 
are likely to have such low [N II] /[O II] ratios, but the KD02 [N II] /[O II] diagnostic can not provide a more specific 
estimate. 

The M91 calibration of R23 makes use of detailed H II region models based on the photoionization code CLOUDY 
(Ferland & Truran 1981). The M91 diagnostic includes the effects of dust and variations in the ionization parameter. We 
have used the analytic expressions for the M91 models given in Kobulnicky, Kennicutt & Pizagno (1999). An initial guess 
is required to determine which branch of the M91 R23 curve to use. We use the [N Il]/[0 II] diagnostic to provide this 
initial abundance estimate. 

The Z94 calibration of R23 is an average of the three independent calibrations given by Edmimds & Pagel (1984); 
Dopita & Evans (1986); McCall, Rybski & Shields (1985), with the uncertainty reflecting the difference among the three 
determinations. A solution for the ionization parameter is not explicitly included in the Z94 calibration. The Z94 
diagnostic was calibrated against H II regions spanning the metallicity range log(0/H) -I- 12> 8.4. As a result, the Z94 
calibration does not reflect the fact that R23 is double- valued with abundance: the use of the Z94 diagnostic assumes that 
all objects have log(0/H) -I- 12> 8.4. We use the [N II]/[0 II] ratio to provide an initial guess of the abundance to ensure 
that the Z94 caHbration is not applied to the objects with log(0/H) -|- 12< 8.4. 

COl gives a number of calibrations depending on the availability of observations of particular spectral lines. Their 
calibrations are based on a combination of stellar population synthesis and photoionization codes with a simple dust 
prescription, and include ratios to account for the ionization parameter. We use the COl "case F" diagnostic which is 
based on the [O III]/H/3 ratio for abundance sensitivity and the [O II] /[O III] 5007 ratio for ionization parameter correction. 
COl recommends using the "case F" diagnostic when the only available emission lines are [O II] [O III], and H/3. For 
[O II]/[0 III]5007< 0.8, the COl diagnostic uses both [O III]/H/3 and [O II]/[0 III]5007, but for [O II]/[0 III]5007> 0.8, 
only the [O III]/H/? ratio is utihzed. Only one of our galaxies has [O II]/[0 III]5007< 0.8, so the COl "case F" diagnostic 
is based on the [O III]/H/? ratio for the majority of our sample. An [O Il]/[0 III]5007 ratio > 0.8 is not unusual: the 
majority of H II regions in van Zee et al. (1998), Kennicutt & Garnett (1996), Walsh & Roy (1997), and Roy & Walsh 
(1997) have [O Il]/[0 III]5007 > 0.8 (Dopita et al. 2000). The COl "case F" diagnostic is potentially problematic for 
our sample because the [O III]/H/3 ratio is relatively insensitive to metallicity (e.g. Dopita et al. 2000). Nevertheless, we 
include the COl "case F" diagnostic because the various COl diagnostics are becoming widely used. 

Figures (6a-d) show the relationship between the metallicity in units of log(0/H)-|-12 and [O II]/Ha (corrected for 
reddening) for the KD02, Z94, M91 and COl abundance diagnostics respectively. The absolute values of the abundances 
vary depending on the diagnostic (Kewley & Dopita 2002). The mean abundances are: log(0/H) -I- 12~ 8.63 (M91), 
~ 8.73 (KD02 [N II]/[0 II]), - 8.60 (COl), and - 8.86 (Z94). Note that the Z94 diagnostic is an overestimate because 
Z94 abundances cannot be calculated for galaxies with log(0/H) -I- 12< 8.4. The KD02 [N II]/[0 II] diagnostic is also an 
upper limit because of the decreasing sensitivity of [N II] /[O II] with smaller abimdances. 

For metallicities log(0/H) + 12> 8.4 (M91, Z94, COl methods) and log(0/H) + 12> 8.5 (KD02 method), we fit a 
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Fig. 6. — Abundance log(0/H) + 12 versus the [O II]/Hq ratio. Abundances are calculated according to: (a) the Kewlcy & Dopita (2002) 
(KD02) [N II]/[0 II] method, (b) the McGaugh (1991) (M91) R23 method, (c) the Zaritsky, Kennicutt & Huchra (1994) (Z94) R23 method, 
and (d) the Chariot & Longhetti (2001) (COl) "case F" method. The dotted lines show the best-fit to the data. For the KD02 diagnostic, 
we fit only data to the right of the dashed line because the KD02 diagnostic becomes less sensitive to abundance below log(0/H) + 12< 8.5. 
For the M91 diagnostic, we fit only data to the rigfit of the dashed line because the M91 diagnostic is double-valued with a local maximum at 
log(0/H) + 12= 8.4. Abundances were not calculated with the Z94 diagnostic for [N II]/[0 II] log(0/H) + 12< 8.4 because Z94 is unreliable 
at these abundances. Errors in the emission-line fluxes are ~ 12%, as described in Kewlcy ct al. (2002). The error in tfie abundance estimates 
is dominated by the inaccuracies of the models used to derive the abundance diagnostics. These errors are ^0.1 dex for the Z94, COl, and 
KD02 diagnostics, and ~ 0.15 dex for the M91 method. The actual error varies depending on the abundance range and the method used, as 
discussed in Kewley & Dopita (2002). 
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least-squares line of best fit to the relationship between [O II] /Ha and abundance (dotted line in Figure 6). This hne has 
the form: 

[OII]/Ha = a * [log(0/H) + 12] + 6 (7) 

where a is the slope and b is the y-intercept. Table 2 gives the slope, y-intercept, and rms for each of the four abundance 

diagnostics. Ideally, all diagnostics should produce the same estimate for the oxygen abimdancc for each galaxy. Unfor- 
tunately, abundance diagnostics are subject to systematic errors resulting from either the modeling, or the data used to 
calibrate the diagnostic (see Kewley & Dopita 2002, for a review). These errors are particularly significant for the R23 
and [O III]/H/3 diagnostics: R23 and [O III]/H/? are double valued with abundance and are strongly influenced by the 
ionization parameter. Because of these errors, the observed relationship between [O II] /Ha and abundance is influenced 
by the shape of the model curves used to calibrate the diagnostics. The shape of the model curves demonstrate the 
theoretical temperature sensitivity of [O II] with increasing abimdance. Because the linear relations are model-dependent, 
it is rather remarkable that for log(0/H) + 12> 8.5, the KD02 [N II]/[0 11], M91 R23, and Z94 R23 diagnostics have 
the same slope and y-intercept to within the errors. Table 2 lists the Spearman Rank coefficients and probabilities. For 
the M91, Z94 and KD02 methods, the correlation coefficient is ~ 0.79 — 0.93, with the probability of obtaining such a 
correlation coefficient by chance < 10""'^^. Because the M91, Z94, and [NII]/[OII] diagnostics arc independent, wc can 
be reasonably confident that the strong correlation observed between [O II] /Ha and abundance is real. Indeed, the R23 
ratio is a valid abundance diagnostic for precisely this reason. For the ionization parameter range of our sample, the 
shape of the R23 curves derives from the temperature sensitivity of the [O II] emission-line compared to H/3. At high 
metallicities, R23 is strongly sensitive to the metallicity because the [O II] and (to a lesser degree) [O III] fluxes drop 
dramatically with the low electron temperatures associated with the increasing abundance. However, at low metallicities 
(log(0/H) -|- 12< 8.4), the electron temperature is high and the [O II] flux increases slowly with abundance. The strong 
relationship between [O II] /Ha and metallicity should also be observed for metallicities derived from non-R23 methods. 
Figure (6a) supports this statement: the [O II] /Ha ratio is strongly correlated with the oxygen abundance derived from 
[N II] /[O II]. The Spearman- Rank correlation coefficient is -0.79 and the probability of obtaining this value by chance is 
negligible (1.75x 10~^^). The slope and y-intercept for the [N II]/[0 II]-derived abundances are within the error range for 
the other three diagnostics. 

The COl [O IIl]/H/3 diagnostic, however, shows a considerably larger scatter, with an rms of 0.24 about the best- 
fit line. KD02 showed that the COl "case a" ([N II] /[S II]) diagnostic also exhibits a larger scatter compared to the 
M91, Z94, or KD02 theoretical methods (including, but not limited to R23). In addition to placing most galaxies at 
abundances log(0/H) -I- 12> 8.5, the COl diagnostic also predicts that 16 galaxies have very low global abundances 
(log(0/H) + 12< 8.0). The COl [O III]/H/3 diagnostic appears to introduce a strong systematic effect in the abundance 
estimates. We will analyze this issue for the COl diagnostic in Section 5. 

4.4. SFR([0 Ilj) and Oxygen Abundance Correction 

In this section, we apply an abundance correction to our [O II] SFR calibration. SFR([0 II]) is normally calibrated 
using an assumed [O II] /Ha ratio. This ratio is not independent of abundance. Wc have shown that the actual [O II] /Ha 
ratio varies considerably for the NFGS sample, and that this variation is strongly correlated with the oxygen abundance. 
For log(0/H) -|- 12> 8.5, the use of any particular [O II] /Ha ratio automatically implies a metallicity which may or may 
not be appropriate for the sample being studied. 

Ideally, one should use the [O II]/Ha ratio for each galaxy to derive an SFR([0 II]) diagnostic. However, if Ha were 
available, it would be used as an SFR diagnostic rather than [O II]. For redshifts z > 0.4, it is theoretically possible to use 
H/3 as a SFR diagnostic through the SFR(Ha) calibration. In practice, H/3 is often contaminated by an unknown amount 
of underlying stellar absorption. In the absence of Ha or a high S/N, high resolution H/3 an oxygen abundance estimate 
can be used as a tracer of the [O Il]/Ha ratio. To obtain the SFR([0 II]), we start with the SFR(Ha) calibration derived 
by K98: 

SFR(Ha)(M0yr"i) = 7.9 x IQ-'^^LiRa) (ergss^^). (8) 
Substituting equation (7) into equation (8) with [O II]/Ha= L(Ha)/L([OII]) gives 



a(log(0/H) + 12) + & 

The SFR([0 II]) spans four orders of magnitude and is therefore particularly sensitive to the values of a, b and log(0/H) -I- 
12. Care should be taken to use a, 6, and log(0/H) + 12 derived from the same abundance diagnostic (Table 2). This process 
assumes (1) that the relationship between [O II]/Ha and metallicity is linear, and (2) that the abundance diagnostic being 
applied is reliable. Both of these assumptions are only valid for metallicities log(0/H) -|- 12 > 8.5 where the [O II] /Ha 
emission decreases with increasing abundance. The [O II] flux is not a strong function of electron temperature at low 
metallicities because the nebular cooling is dominated by hydrogen free-free emission. 

Figure 7 shows the relationship between the K98 Ha SFR and the SFR([0 II]) derived from our new calibration 
(equation 9) for each of the abundance diagnostics. In each plot, a dotted line indicates the best flt to the data. Table 3 
gives the slope, y-intercept, and rms for each fit. For comparison, Table 3 also lists the slope, y-intercept, and rms for 



12 



the K98 SFR([0 II]) and SFR(Ha) plot in Figure 2. After correction for oxygen abundance, in all four cases the line of 
best fit to the data has a slope of ^ 1 and a y-intercept of ~ within the errors, indicating that the abundance correction 
does not introduce a systematic offset. For the KD02, M91 and Z94 diagnostics, the rms scatter decreases significantly 
after correction for oxygen abundance (0.03-0.05 versus 0.08-0.11). 

Cardiel et al. (2003) also observed a decreased scatter after metallicity correction in a small sample of 7 galaxies with 
redshifts of z ~ 0.4 and z ~ 0.8. Cardiel ct al. applied a metallicity correction to [O II] based on R23 and found 
excellent agreement between SFR([0 II]) and SFR(Ha). This result gives us confidence that our abundance-corrected 
SFR([0 II]) calibration will be applicable to more distant samples than the NFGS. Indeed, we derive our SFR([0 II]) 
calibration only from the strong [O II] /Ha-metallicity correlation. In theory, this correlation is a result of the temperature 
sensitivity of [O 11] relative to Ha and, therefore, should not be sensitive to redshift. In practice, however, the situation 
is more complicated. The abundance diagnostics are based on theoretical models calibrated against nearby H II regions 
or galaxies. It is unclear whether the model assumptions apply at high-2;. Model assumptions which may differ at high-2; 
include (but are not limited to) the gas geometry, dust geometry, density, and the initial mass function. 




Fig. 7. — Bottom panels: SFR(Ha) versus SFR([0 II]) for the four abundance diagnostic methods. The SFR([0 II]) ratio is corrected for 
reddening and abundance according to equation (10). Abundances arc calculated using: (a) the Kewley & Dopita (2002, ; KD02) [N II]/[0 II] 
method, (b) Zaritsky, Kennicutt & Huchra (1994, ; Z94), (c) McGaugh (1991, ; M91), and (d) the Chariot & Longlietti (2001, ; COl) "case 
F" method. Top panels: The K98 SFR(Ha) versus the logarithm of the ratio of SFR([0 II]) and SFR(Ha) from the bottom panels. In each 
panel, the dotted line shows the least squares fit to the data. Estimated errors are ~ 30% for SFR(Ha) and ~ 35% for SFR([0 II]), as in 
Kewley et al. (2002). 
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The large scatter (Figure 7d) for the COl "case F" ([O III]/H/3) abundance diagnostic propagates into the SFR([0 II]) 
cahbration based on the COl constants and abundance estimate. Wo therefore do not recommend the use of COl case F 
to derive an abundance-corrected [O II] star formation rate if large scatter is a concern. As we have seen, the M91, Z94 
and KD02 abundance diagnostic methods (and associated a and b constants) give almost identical relations (within the 
errors) between SFR(Ha) and SFR([0 II]) with a very small scatter (0.03-0.05 dcx). The fact that Z94, M91, and KD02 
([N II] /[O II]) arc independent of one another and still produce identical relations (within the errors) supports the use of 
these to correct the NFGS [O II] SFRs for oxygen abundances between log(0/H) + 12~ 8.5 - 9.2. 

The drawback to using R23 diagnostics is that they are double-valued with abundance. The M91 diagnostic requires 
an initial guess of the oxygen abundance to determine which branch of the R23 curve to use. The Z94 calibration is only 
valid for the upper mctallicity branch (log(0/H) + 12> 8.4). Unfortunately, the [O II], [O III], and H/3 lines alone are not 
sufficient to determine which branch of the R23 curve to use. For the NFGS sample, we use the [N II] /[O II] line ratio 
to resolve this problem. In local galaxies, the luminosity-metallicity (L-Z) correlation may help to break the degeneracy. 
For example, objects more luminous than Mb — —18 generally have metallicities greater than log(0/H) -|- 12~ 8.4 (e.g., 
Z94) and therefore probably lie on the upper R23 branch. Figure 8 supports this conclusion. In Figure 8, we compare 
the absolute magnitude Mb for the NFGS galaxies with the abimdanccs derived with the KD02 [N II] /[O II] and the 
M91 R23 methods. Even though the KD02 [N II]/[0 II] method is less sensitive to abundance for log(0/H) + 12< 8.5, 
the KD02 [N II] /[O II] method shows a strong correlation between abundance and Mb. For abundances estimated using 
KD02, all eight galaxies with log(0/H) + 12< 8.4 have Mb < -18. For abundances calculated using M91, 13/15 (87%) of 
the galaxies with log(0/H) + 12< 8.4 have Mb < —18. For Mb < —18, therefore, Mb is a useful discriminator between 
the two R23 branches in nearby galaxies and provides a crude estimate of the abundance in the absence of alternative 
methods. The error in the abundance is likely to be ±0.2 in units of log(0/H)-|-12. At lower luminosities (Mb > —18), 
the Me-metallicity relation provides, at most, an upper limit. 

It is not clear whether the same MB-metallicity relationship applies for galaxies at higher redshifts. The few studies of 
the luminosity-metallicity (L-Z) relation at larger redshifts appear to produce conflicting results. CaroUo & Lilly (2001) 
analysed a sample of 13 star forming galaxies between 0.5 < z < 1 and find no significant evolution in the L-Z relation 
out to 2; = 1. Lin et al. (1999) examined >2000 late-type CN0C2 (Canadian Network for Observational Cosmology 
Field Galaxy Redshift Survey) galaxies and found no significant luminosity evolution between 0.12 < z < 0.55. However, 
results from the DEEP Groth Strip Survey suggest that the L-Z relation does evolve from the local relation between 
z = to z = 1 (Kobulnicky et al. 2003). At larger redshifts z > 2, the L-Z relation appears to be significantly different 
from the local relation (Kobulnicky & Koo 2000; Pettini et al. 2001). Therefore, although potentially useful, the local 
MB-metallicity relationship should not be applied blindly to non-local samples. 

To conclude, the Z94 abundance estimates agree well with those obtained using the KD02 [N II] /[O II] diagnostic. The 
[N II] /[O II] ratio is very sensitive to mctallicity and is almost independent of ionization parameter (KD02). Therefore, if 
an initial guess of the abundance gives log(0/H) -|- 12> 8.4, and the ionization parameter is likely to cover a small range 
(similar to the NFGS or local H II regions), we recommend using the Z94 R23 diagnostic. Using the slope and y-intercept 
appropriate for Z94 (Table 2) gives: 

QPpanTTl -1^ 7.9 X 10-4^ L([OII]) (ergs s^i) 

SFR([OII],Z)(Meyr ) = ± o.25) . [log(0/H) + 12] + (16.73 ± 2.23) ' ^'^^ 

where log(0/H) -|- 12 comes from 



log(0/H) -I- 12 = 9.265 - 0.33 R23 - 0.202 R^g - 0.207 R^g - 0.333 R^g (Z94) (11) 

and R23 = log ([Oil] A3727 -I- [OIII] A5007)/H/3. 

Given the difficulty in estimating abundances with limited data, our SFR([0 II], Z) calibration should be useful for 
deriving an SFR([0 II]) calibration for samples which have a different mean abundance from the NFGS. Equation (4) is 
based on the mean intrinsic [O II] /Ha of the NFGS sample. Any SFR([0 II]) calibration that is derived from SFR(Hq:) and 
an [O II] /Ha ratio automatically includes an assumption about the average abundance. As we have seen, the remaining 
cause of discrepancies in [O II] /Ha from galaxy to galaxy in the NFGS is abundance. If the mean abundance of a sample 
is not the same as for the NFGS (using the same abundance diagnostic), then equation (10) can be used to calculate a 
new SFR([0 II]) calibration based on the mean or assumed abundance for the new sample. This approach could be useful 
in cases where individual galaxy abundances are not available, but an estimate of the sample mean abundance can be 
made. Such estimates could be based on known abundances for similar galaxies, or could be calculated using a subsample 
of galaxies for which abundance measurements are available (e.g., Lilly, CaroUo, & Stockton 2003). A similar process can 
be utilized for deriving a mean sample extinction estimate. 

If no abundance estimate can be made, we recommend using equation (4) derived in Section 4.1. This SFR indicator 
is most useful for large samples because, provided the mean abimdance is similar to that observed in the NFGS, the 
mean SFR([0 II]) should approximate the mean SFR(Ha), thus reducing the scatter. Note that our equations (4) and 
(10) assume that the [O II]/Ha ratio does not depend significantly on the ionization parameter. Investigations using 
SFR([0 II], Z) should also include a calculation of [O III]/[0 II] to measure the dominant ionization state of oxygen. If 
[O in]/[0 II] covers a wide range, and the oxygen abundance covers a relatively small range, then equation (4) would be 
a more appropriate SFR([0 II]) calibration to use. 
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5. A THEORETICAL CALIBRATION OF SFR([0 II]) AND ABUNDANCE 

In this section, we utilize theoretical models to further investigate the relationship between the [O II] /Ha ratio, abun- 
dance, and the ionization state of the gas. Wc use the stellar population synthesis models Pcgasc (Fioc & Rocca- 
Volmerange 1997) and StarburstQQ (Leitherer et al. 1999) to provide the ionizing stellar radiation field for the photoion- 
ization code, Mappings III (eg., Sutherland & Dopita 1993; Groves et al. 2003). Mappings III self-consistently calculates 
radiative transfer through gas in the presence of dust. Our models, described in Kewley et al. (2001b); Dopita et al. 
(2000), have been successfully applied to H II regions (Kewley & Dopita 2002; Dopita et al. 2000) and nearby starburst 
galaxies (Kewley et al. 2001b; Calzetti et al. 2003). Wc use the instantaneous burst models with an ionization parameter 
range of g = 5 x 10^ — 8 x lO^cm/s. The models cover metallicities of 0.05, 0.1, 0.2, 0.5, 1.0, 1.5, 2.0, and 3.0xsolar, 
where solar metallicity is defined in Anders & Grevesse (1989). The corresponding metallicities in log(0/H) + 12 are 
7.6, 7.9, 8.2, 8.6, 8.9, 9.1, 9.2, 9.4. Note that for the currently favored value of solar abundance (log(0/H) + 12 - 8.69; 
AUende Prieto, Lambert, & Asplund 2001), the model metallicities become 0.09, 0.2, 0.4, 0.9, 1.7, 2.6, 3.5, 5.2xsolar. The 
metallicities correspond to specific stellar tracks used in the population synthesis models and to the nebular abundance of 
the photoionization models. Typical metallicities for H II regions range between 8.2 < log(0/H) + 12 < 9.2 (e.g., Kewley 
& Dopita 2002). 

Figure 9 shows the theoretical relationship between [O 11] /Ha and metallicity. We conclude that: (a) both ionization 
parameter and metallicity affect the [O II] /Ha ratio, and (b) for a particular sample, the relative importance of the 
ionization parameter compared to metallicity is governed by the range in abundances and ionization parameters spanned 
by the sample. The [O II] /Ha ratio is relatively insensitive to metallicity for 8.3 <log(0/H) + 12< 8.6. However, for 
log(0/H) + 12> 8.6, [O II]/Ha becomes a strong function of metallicity, particularly for low ionization parameters. This 
behaviour reflects the temperature sensitivity of the [O II] emission-line. As wc have discussed, for temperatures typical 
of star-forming regions (10000-20000 K), the excitation energy between the two upper D levels for [O II] and the lower 
S level is of the order of the thermal electron energy kT. The [O II] doublet is therefore closely linked to the electron 
temperature. At low metallicities, the electron temperature is high, and [O II] emission increases with metallicity. In 
this regime, the thermal cooling is dominated by hydrogen free-free emission. However, when the metallicity increases, 
the number of coolants in the nebula rises, thus lowering the electron temperature. The [O II] emission therefore drops 
rapidly with increasing metallicity. 

Figure 9 predicts that samples covering a small range of metallicities will not show a correlation between [O II] /Ha and 
abundance because of the range of possible ionization parameters. Samples spanning metallicities 8.3 <log(0/H)-|-12< 8.6 
will also not show a strong relationship between [O II]/Ha and log(0/H) -|- 12 because at these metallicities [O II]/Ha 
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Fig. 8. — comparison between the metallicity in units of log(0/H) -|- 12 and the blue absolute magnitude. Abundances are calculated using: 
(a) the Kewley & Dopita (2002, ; KD02) [N II]/[0 II] method, and (b) the McGaugh (1991, ; M91) R23 method. The error in the blue 
absolute magnitude is 0.013 mag on average (Jansen et al. 2000a). The error in the abundance estimates is < 0.1 dex for the Z94, COl, and 
KD02 diagnostics, and ~ 0.15 dex for the M91 method. The actual error varies depending on the abundance range and the method used, as 
discussed in Kewley & Dopita (2002). 
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is a weak function of log(0/H) + 12 and [O II]/Hq; is more strongly affected by the ionization parameter. Only samples 
with a normal range of ionization parameters = 1 x 10^ — 8 x lO^cm/s; Dopita et al. 2000) and covering a large range 
of metallicities exceeding log(0/H) + 12'^ 8.5 will exhibit a correlation between [O Ilj/Ha and metallicity. 

Note that our models do not require a specific method for calculating metallicities from the data. If the metallicities are 
calculated using some reliable abundance diagnostic, our models predict that galaxies with typical ionization parameters 
and metallicities lie along the curves in Figure 9. The curve for each ionization parameter can be characterized by a third 
order polynomial: 



[Oil] 
Ha 



bx 



(12) 



where x =log(0/H) + 12 and the coefficients a, 6, c, d arc displayed in Table 4. 

In Figure (10), we compare the theoretical models with the abundances derived using the different diagnostics. In 
Figure (10a) we show the abundances derived using [N II] /[O II]. The ionization parameters are typically 2 x 10^ — 4 x 
10'' cm/s , and the data follow a similar trajectory to the models. As we have discussed, the oxygen abundances derived 
using the [N II]/[0 II] diagnostic show a larger scatter, particularly for oxygen abundances log(0/H) + 12< 8.5. The 
models show that the [N II] /[O II] and [O II] /Ha ratios become less sensitive to abundance as metallicity decreases, 
increasing the scatter. The Z94 method produces abundances similar to the [N II] /[O II] method. The data follow a 
similar trajectory to the models and the ionization parameters are between typically 2 x 10^ — 4 x 10^ cm/s. 

The M91 method (Figure 10c) shows a systematic offset compared to the Z94 and KD02 methods and to the models. This 
offset has been observed previously (Kewley & Dopita 2002) and is probably a result of the different stellar atmospheres 
and stellar models used to derive the ionizing radiation fields. The difference in abundances estimated for the M91 and 
KD02 method is ~ 0.1 — 0.2 in log(0/H) + 12. This variation is within the errors associated with the diagnostics (0.15 
dex for M91 and 0.1 for KD02). We note that the R23 diagnostics may be minimizing the scatter by making a hidden 
assumption about the ionization parameter. The R23 ratio is sensitive to the ionization parameter and the ionization 
parameter diagnostic [O III]/[0 II] is sensitive to abundance. If the ionization parameter correction is not made iteratively, 
then a calibration may, in effect, favor a particular ionization parameter or range of ionization parameters. This effect 
may contribute to the difference between the M91 and Z94 diagnostics. The Z94 abundance estimates agree well with the 
ionization-parameter independent diagnostic [N II] /[O II]. 

The COl "case F" diagnostic produces some [O II]/Ha-abundance combinations which can not be produced using our 
stellar population synthesis+photoionization models, even with the 0.1 dex error estimates. As we discussed earlier. 



O 1.0 




8 6 8 8 

log(0/H) + 12 (model) 

Fig. 9. — Theoretical models of [O II] /Ha as a function of the model gas-phase oxygen abundance compared to the NFGS galaxies. The 
colored lines correspond to different ionization parameters, shown in the legend in cm/s. The model grid errors are < 0.1 dex in log(0/H) + 12, 
and < 23% in [O II] /H«. 
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Fig. 10. — Theoretical models of [O IIJ/Ho as a function of the model gas-phase oxygen abundance compared to the NFGS galaxies. 
The model oxygen abundance is defined by the metallicity of the stellar tracks in the stellar population synthesis models, and is the same 
motallicity used for the surrounding nebular gas. The NFGS abundances were calculated using: (a) the Kewley & Dopita (2002, ; KD02) 
[N II]/[0 II] method, (b) Zaritsky, Kennicutt & Huchra (1994, ; Z94), (c) McGaugh (1991, ; M91), and (d) the Chariot & Longhetti (2001, ; 
COl) "case F" method. 
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the COl "Case F" diagnostic is based on [O III]/H/? for the majority of galaxies in our sample. Figure 11 shows the 
theoretical relationship between [O III]/H/3 and abundance. The [O III]/H/3 ratio is much more sensitive to the ionization 
parameter than metallicity for all but the highest metallicities. In addition, [O III]/H/3 is double valued with abundance. 
Any particular [O III]/H/3 ratio could correspond to a range of abundance/ionization parameter combinations, and the 
possibility of obtaining an incorrect abundance estimate is high. 

We can conclude from Figures lOa-c that the [O II] /Ha ratio depends on abundance for log(0/H) + 12> 8.5, and that 
a linear correction for abundance is theoretically plausible for samples with metallicities in this range, providing the data 
do not span a larger range in ionization parameters than is observed in local H II regions. For the KD02, Z94, and M91 
abundances, the NFGS data follow a trajectory with a similar slope to the models for log(0/H) + 12> 8.5. The model 
trajectory for each ionization parameter can be used to derive a theoretical [O II] SFR calibration. We begin, once again, 
with the K98 SFR(Hq!) calibration: 



SFR(Ha) (Mgyr-i) = 7.9 x 10~^^L(Ra) (ergss-^) 
_ 7.9 X 10-42 L([OII]) (ergs s-i) 
^ [Oil] /Ha 

Substituting equation (12) into equation (13) for [O II]/Ha then yields: 



(13) 



sFRdoiu.z), = "°'"°:r^™^T"'' (») 

a + ox + cx'^ + dx-^ 

where the constants a,b,c,d are given in Table 4 for each ionization parameter and x =log(0/H) + 12. The subscript 
t indicates that the calibration is based on theoretical models. The majority of the NFGS have ionization parameters 
between q = 2 x 10^ — 4 x lO^cm/s, according to the KD02 [N II]/[0 II] and Z94 R23 diagnostics. Interpolating between 
the q ~ 2 X 10^ and 4xl0^cm/s curves gives a curve with an approximate ionization parameter oi q — 3 x lO^cm/s: 

SFR([OII],Z), (M,yr-) = 7.9 x IQ-^^ L([OII]) (ergs s-^) 

^' ^* ^ "■^ ' -1857.24+ 612.693a; - 67.0264x2 + 2.43209x3 ^ ' 



1.01 ' ' ' ' — \ ' ' ' ' — \ ' ' ' ' — \ ' ' ' — r 




LOG (0/H) + 12 

Fig. 11. — Theoretical models of [O III]/H/3 as a function of the model gas-phase oxygen abundance. The model oxygen abundance is 
defined by the metallicity of the stellar tracks in the stellar population synthesis models, and is the same metallicity used for the surrounding 
nebular gas. The [O III]/H/3 ratio is almost independent of metallicity for all but the highest metallicities. The grid errors are estimated to 
be < O.ldex in log(0/H) + 12, and up to 23% in [O III]/H/3. 
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This curve provides a useful theoretical description of the behavior of [O II] /Ha with mctallicity for the NFGS sample. 
We emphasize that the metallicitics of the models correspond to the mctallicities in the stellar tracks and the modeled 
nebulae, and are independent of the method used to derive log(0/H) + 12. Therefore, any method can be used to derive 
log(0/H) + 12, as long as the method is reliable over the expected abundance range of the sample. The Z94 R23 diagnostic 
can easily be used if the abundances exceed log(0/H) + 12> 8.4. 

In Figure 12, we compare the Ha and [O II] SFRs with SFR([0 II], Z) calculated according to our theoretical models 
(equation 15) with abundances estimated by either (a) the KD02 [N II] /[O II] method or (b) the Z94 R23 method. Table 3 
contains the slope, intercept, and scatter. The slope and y-intercept are close to unity and zero respectively, for both 
abundance diagnostics. The scatter is 0.05 dex compared to 0.08 dex using our SFR([0 II]) calibration without correction 
for abundance (equation 4; Figure 4). Clearly our theoretical SFR([0 II],Z)i calibration is successful in reducing the scatter 
observed in the SFR([0 II]) estimates. This diagnostic will be most useful for deriving a new SFR([0 II]) calibration for 
samples which have a different mean abundance from the NFGS. 

6. THE APPLICATION OF SFR([0 II]) TO HIGH Z GALAXIES 

6.1. Reddening Determination 

Recently, many investigations have used [O II] to constrain the cosmic star formation history for redshifts 0.4 < 2; < 1.6 
(e.g.. Hammer et al. 1997; Hogg et al. 1998; Rosa- Gonzalez, Terlevich, & Terlevich 2002; Hippelein et al. 2003). At these 
redshifts. Ha is usually unavailable and correction for reddening using the methods outlined above is thus impossible. 
Without the Balmcr decrement, many investigators apply an "average" or "recommended" mean attenuation of Ay ~ 
1 mag prior to the calculation of either SFR([0 II]) or abundance. Assuming Ry = Ay/E(i3 — V) = 3.1, Ay = 1 
corresponds to E(B — V) ~ 0.3. Figure 13 shows the distribution of reddening traced by E(B-V) for our sample. The 
mean E(B-V) for the galaxies in our sample (after correction for Galactic extinction) is 0.26 ± 0.02, consistent with the 
common choice of E(_B — y) ~ 0.3. 

If we apply an E(_B — V) — 0.3 to R23 and L([0 II]) and use equations (lO)-(ll) or equation (15) to derive the SFR, 
the slope is a ~ 0.77 ± 0.03 (Figure 14a,b). Thus, with a single E(j3 — the SFR([0 II]) is a systematic underestimate 
at high SFRs and a systematic overestimate at low SFRs. This effect would be observed if the galaxies at the highest 
SFRs are more highly extinctcd than galaxies with lower SFRs. Wang & Hcckman (1996) showed that the reddening 
(measured using the Ha/H/3 ratio) correlates with FIR luminosity for a sample of nearby disk galaxies. A similar effect 
appears in a sample of nuclear starburst and blue compact galaxies by Calzetti ct al. (1995). We also know from Kewley 
et al. (2002) that the SFR(FIR) agrees to within 10% on average with SFR(Ha), and we have shown here that there 
exists a 1:1 relationship between SFR(Ha) and SFR([0 II]) after reddening and abundance correction. It is therefore not 
surprising that we observe increasing reddening with SFR([0 II]). 




Fig. 12. — Comparison between the SFRs derived from Ha and [O II] with SFR([0 II], Z) calculated using our theoretical grids (equation 15) 
for an eflfective ionization parameter q ~ 3x lO'^cm/s. Abundances were calculated using: (a) the Kewley & Dopita (2002, ; KD02) [N II]/[0 II] 
method, and (b) the Zaritsky, Kennicutt & Huchra (1994, ; Z94) R23 method. The error in SFR(Ha) is ~ 30% and the error in SFR([0 II],Z) 
is ~ 35%. 
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Fig. 13. — The reddening distribution traced by E(B-V) for the NFGS. All fluxes have been corrected for Galactic extinction. The vertical 
line at the top shows the mean of the distribution: 0.26 ± 0.02. 




Fig. 14. Bottom panel: Comparison between the SFR(Ha) and SFR([0 II]) with the average reddening of Av = 1 applied to every 
galaxy. The SFR([0 II] ratio is corrected for reddening and abundance according to (a) equation (10) and (b) equation (15) Abundances are 
calculated using the Z94 R23 method, which is comparable to the [N II] /[O II] method. The solid line has a slope of one; the dotted line is 
the least squares best fit to the data. The estimated error in SFR(Ha) is ~ 30% and the error in SFR([0 II], Z) is ~ 35%. The SFR([0 11], Z) 
error does not include the systematic error introduced by using Av = 1 rather than the reddening derived from the Balmer decrement. Top 
panel: The K98 SFR(Ha) versus the logarithm of the ratio of SFR([0 II]) and SFR(Ha) from the bottom panel. 
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Fig. 15.— Comparison between (a) the L([0 II]) and E(B-V), and (b) the L([N II]) and E(B-V). L([0 II]) and L([N II]) are corrected for 
reddening using the Balmer decrement. The solid horizontal line indicates the detection limit on E(B-V) for our data. Galaxies lying on or 
below this line have upper limits for E(B-V). The dashed line is the least squares fit to the data if the upper limits are fixed at 0.02. The 
dot-dashed line is the least squares fit if the upper limits are fixed at 0.00. Estimated errors in the line luminosities are ~ 16% while the 
errors in E(B-V) are ±0.04. 



Figure 15 sliows a strong correlation of L([0 II]) with E(B-V) for our sample with a large scatter (rms=0.13 dex). We 
corrected L([0 II]) for reddening using the Balmer decrement. The Spearman Rank coefficient is 0.73; the probability 
of obtaining this coefficient by chance is 8 x 10^^''. We made two fits to the data to estimate the impact of the upper 
limits: one with the upper limits on E(B-V) set at their maximum value of 0.02 (dashed line in Figure 15), and the second 
line with the upper limits set at zero (dot-dashed line). The fits are almost identical: the slope and y- intercepts agree to 
within 1%, well within the errors. The best fit (including E(B-V) upper limits as either 0.02 or 0.00) is: 

E{B -V) = (0.174 ± 0.035) log[L([OII])i] - (6.84 ± 1.44) (16) 

where the subscript i indicates that L([0 II]) (in ergs/s) has been corrected for reddening. The intrinsic [N II] luminosity 
is also correlated with E(B-V) (Figure 15). The [N II] emission-line requires less correction for reddening than [O II] and 
provides an independent measure of the relationship between the emission-line luminosity and reddening for the NFGS 
sample. 

This relationship can be understood physically: the galaxies with the highest rates of star formation are likely to also 
produce larger quantities of dust. Most of the dust in galaxies is probably produced by carbon or oxygen-rich stars on the 
asymptotic giant branch (see Mathis 1990, for a review). Supernovae may also be important because they insert heavy 
elements into the surrounding interstellar medium. If we assume that the initial mass function is similar for the galaxies 
in our sample, then higher star formation rates enable more carbon and oxygen-rich stars to reach the asymptotic giant 
branch, producing larger quantities of dust. This scenario implies that most of the NFGS galaxies must have been forming 
stars at least 1.5-2 Gyr ago because it takes approximately this long for the low- and intermediate- mass stars in a typical 
stellar population to evolve to the asymptotic giant branch (e.g., IVIouhcine & Lancon 2002). 

The scatter around the fit in Figure 15 is large (0.13 dex). The amount of dust obscuring the observed optical emission 
from the nebular gas varies from galaxy to galaxy. In addition, because we use global spectra, we observe the sum of the 
emission from the brightest H II regions in each galaxy. Geometry, dust composition, and stellar properties are all likely 
to have an impact on the observed optical emission from the brightest H II regions. 

Because the reddening is correlated with the intrinsic [O II] hmiinosity, equation (16) provides a very crude estimate 
of the reddening for the NFGS. The relationship between E(B-V) and the [O II] luminosity could be different for other 
samples and tests are required to determine whether equation (16) may be applied to non-NFGS galaxies. In addition, 
for galaxies at high redshift E(B-V) may not be a reliable indicator of the reddening if the dust does not conform to 
a foreground screen geometry (Witt & Gordon 2000). However, there is some evidence that a reddening-luminosity 
relationship exists at high redshifts, at least for rapidly star-forming galaxic;s. Adelberger & Steidel (2000) found that the 
sum of the bolometric dust luminosity and the I6OOA luminosity (Lboi.dust + Lieoo) is correlated with the ratio of these 
two quantities (Lboi,dust/Li6oo) for high- 2; and low- 2; galaxies alike. The sum Lboi,dust + Lieoo provides a crude estimate of 
the star formation rate; the ratio Lboi,dust/Li600 is a rough tracer of the dust obscuration. If such a reddening- luminosity 
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relationship holds for non-NFGS samples, then it might be feasible to use equation (16) to derive a very rough reddening 
estimate at distances where Ha is redshifted out of the observable wavelength range. 

The intrinsic [O II] luminosity L([OII])i is related to the observed [O II] luminosity L([OII])o using the standard equation 
(e.g., Calzetti et al. 2000): 

L([OII])i(ergss-i) = L([OII])o(ergs s-^) x io0-4fe[onj e(b-v) ^^j-^ 

where fc[oii] = 4.771 using the COM reddening law. If we substitute equation (16) into equation (17) for E(B-V), we 
obtain: 

L([OII])i = 3.11 X 10-20l([OII])^-^^^ (18) 

where the intrinsic and observed luminosities are in units of ergs/s. The estimated intrinsic luminosity from equation (18) 
can now be used in equation (10) along with R23 to calculate the SFR([0 II], Z) for the NFGS. We emphasize that this 
relation may be different for other samples, and should not be applied blindly to other galaxies. 

Figure 16 compares SFR([0 II] ,Z) with SFR(Ha) for the NFGS, where SFR([0 II], Z) is effectively corrected for reddening 
and abundance in the following manner: 

1. We estimate the intrinsic L([OII])i luminosity from equation (18). 

2. We estimate E(B-V) from the intrinsic L([OII])i luminosity with equation (16). 

3. We use the E(B-V) estimate to correct the R23 ratio for reddening. 

4. We use the reddening-corrected R23 ratio in the Z94 diagnostic (equation 11) to derive the abundance. 

5. We use the abundance and intrinsic L([OII])i luminosity in equation (10) to derive a reddening and abundance 
corrected SFR([0 II], Z) estimate. 

6. We use the abundance and intrinsic L([OII])i luminosity in equation (15) to derive a reddening and abundance 
corrected SFR([0 II],Z)t estimate from our theoretical SFR([0 II],Z)t calibration. 




Fig. 16. — Bottom panel: Comparison between the SFR(Ha) and SFR{[0 II], Z). SFR([0 II], Z) is corrected for reddening and abundance 
using the E(B-V) derived from the correlation of E(B-V) with [O 11] luminosity (equation 16). In panel (a), the abundance was calculated 
using equation (10) (best-fit to the [O II]/Ha and log(0/H) + 12curve). In panel (b), the abundance was calculated using equation (15) 
(derived from our theoretical models) . The dotted line shows the least squares fit to the data. The error in SFR(Ha) is ~ 30% and the error 
in SFR([0 II], Z) is ~ 35%. The SFR([0 II], Z) error does not include the error introduced by using the [O II] luminosity rather than the 
reddening derived from the Balmer decrement. Top panel: The SFR(Ho;) versus the logarithm of the ratio of SFR([0 II]) and SFR(Ha) from 
the bottom panel. 
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Note that we correct the comparison star formation rate SFR(Hq:) for reddening using the true E(B-V) obtained from 
the Balmer decrement. The best-fit hues in Figure 16 have a slope a = 1.07 ± 0.03 and a y-intercept close to zero 
(6 = —0.06 ± 0.03). Clearly, if the best-fit slope is not close to 1 (as is the case when an average A(v)=l is assumed), 
the difference between the mean estimated SFR and the true SFR increases with increasing star formation rate. Thus 

the bias introduced can easily appear as a function of rcdshift because we observe intrinsically more luminous galaxies 
at larger z. The scatter in Figure 16 is 0.21 dcx but this scatter is less important than the slope for studies of the star 
formation history. 

Figure 17 shows a comparison of the SFR(Hq) with the SFR([0 II]) where SFR([0 II]) is efi^ectively corrected for 
reddening but only partially for oxygen abundance through the correlation of E(B-V) with luminosity. We calculated the 

intrinsic L([011])i luminosity using equation (18), and derived the SFR([0 II]) using equation (4). The scatter is larger 
than if we apply an abundance correction (~ 0.26 dex vs. 0.21 dex), but again, with large samples this increased scatter 
should not be a problem. 

As we have stressed, our equations 4 and 10 were derived from an unbiased local sample. The excitation and abundance 
properties of galaxies in other samples and of galaxies at higher redshifts may not be the same as those observed in the 
local galaxy population. It may be possible to correct L([OII])o for reddening using a reddening- huninosity relation such 
as equation 18, but the application of equation 18 (or other such relation) awaits further testing. The lack of samples with 
both integrated spectra and Balmer decrement measurements (or other reddening indicator) make testing equation 18 
difficult at present. We will test the reddening-luminosity relation for a large, objectively selected sample of galaxy pairs 




Fig. 17. — Bottom panel; Comparison between the SFR(Ho) and SFR([0 II]). SFR([0 II]) is corrected for reddening using the E(B-V) 
derived from the correlation of E(B-V) with [O II] luminosity (equation 16). SFR([0 II]) is not explicitly corrected for abundance. SFR([0 II]) 
is only corrected for abundance through the correlation of E(B-V) with L([0 II]). The dotted line shows the least squares fit to the data. The 
error in SFR(Ha) is ~ 30% and the error in SFR([0 II],Z) is ~ 35%. The SFR([0 II],Z) error does not include the error introduced by using 
the [O II] luminosity rather than the reddening derived from the Balmer decrement. Top panel: The SFR(Ha) versus the logarithm of the 
ratio of SFR([0 II]) and SFR(Ha) from the bottom panel. 
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in a future paper (Kcwley et al. in prep), once aperture effects liave been analysed. For the reasons outlined above, we do 
not apply equation f8 to high-z galaxies. Instead, we test whether our SFR([OII]) relations can be applied to galaxies at 
higher redshifts with the standard assumption of Ay = 1. Evidence for a mean Ay ~ 1 for galaxies with redshifts z < 0.3 
is found by Tresse et al. (2002), Sullivan et al. (2000), Liang et al. (2003) using the Balmer decrement. However Flores et 
al. (1999) concludes that the global opacity of the {z < 1) universe is between Ay = 0.5 — 0.85 using radio, IR, UV, and 
optical photometry. 

6.2. Testing SFR([OII]) for galaxies 0.5 < z < 1.6 

To test our SFR([OII]) relations on galaxies at larger-z, we use two samples: the sample of Hicks et al. (2002) (hereafter 
H02) and the sample of Tresse et al. (2002) (hereafter T02). H02 obtained rest-frame blue spectra for 14 emission-line 
galaxies representative of the population at 0.8 < z < 1.6 in the 1999 NICMOS parallel grism Ha survey. We use the seven 
objects in H02 with measured [O II] and Hq luminosities. T02 obtained Ha measurements for 30 galaxies in the Canada- 
France Redshift Survey with redshifts 0.5 < 2 < 1.1. Both H02 and T02 report [O II] /Ha ratios that differ significantly 
from those in the NFGS, due either to much larger reddening values or intrinsically lower line ratios. Therefore, these 
samples provide a critical test of the applicability of our [O II] -SFR calibration to galaxies at higher redshifts. 

Wc recalculated the [O II] and Ha luminosities for the standard cosmology {h = 0.72; flm = 0.29) and computed 
the intrinsic [O II] and Ha luminosities with the CCM reddening curve, assuming Ay — 1- Table 5 lists the relevant 
derived quantities. The SFR(Ha) and SFR([0 II]) are from equations 1 and 4. SFR([0 II]) is significantly lower than 
SFR(Ha) (Figure 18a). This discrepancy may result from one or a combination of (1) poor (< 5a) S/N at [O II], (2), sky 
contamination at [O II], (3) a stronger dependence of [O II] /Ha on the ionization parameter than observed in the NFGS, 
(4) reddening more than the assumed Ay = 1, (5) reddening may not be a simple forgroimd screen for the H02 and T02 
galaxies, (6) the H02 and T02 samples may have a different mean oxygen abundance than the NFGS. We discuss these 
possibilities individually below. 

6.2.1. Poor S/N 

Poor S/N at [O II] may affect the H02 sample: only two H02 galaxies have > 5a [O II] detections. However, the 
majority (28/30) of the T02 sample have > 5a detections, ruling out poor S/N as a reason for the SFR discrepancy for 
the majority of galaxies in Figure (18a). 

6.2.2. Sky ContaminaMon 

Sky contamination at [O II] affects the majority (5/7 galaxies) of the H02 sample, but the T02 sample was specifically 
selected to avoid sky contamination in the optical and near-IR. Inspection of the T02 spectra reveals no evidence for sky 
contamination. 

6.2.3. Ionization Parameter 

Ionization parameter estimates can not be made with the T02 and H02 samples because [O III] A5007 emission-line 
fluxes are unavailable. However, [O II] and [O III] fluxes arc available for the 0.47 < z < 0.92 sample by Lilly, Carollo, & 
Stockton (2002, 2003) (hereafter LCS). The mean [O III]/[0 II] ratio for the LCS sample is 0.37 ± 0.04, corresponding to 
a mean ionization parameter of ~ 3.5 x lO'' cm/s using the [O III]/[0 II] ionization parameter calibration of Kewley & 
Dopita (2002). If the mean ionization of the H02 and T02 sample is similar to that observed by LCS, then the ionization 
parameter cannot account for the SFR discrepancy in Figure (18a). 

6.2.4. Reddening 

If the SFR discrepancy is a result of reddening, then the average extinction required to bring the SFRs into agreement 
is Ay ~ 1.6. The variation in optical extinction as a function of redshift is unknown, however the multiwavelength study 
by Flores et al. (1999) suggests a mean extinction in < 2: < 1 galaxies of Ay = 0.5 — 0.85, significantly lower than 
Av ^ 1.6. 

6.2.5. Foreground Screen Assumption 

An alternative explanation for the difference between the [O II] and Ha SFRs in Figure (18a) is that the dust geometry 
may not be a simple screen. Witt & Gordon (2000) show that E{B — V) saturates at around 0.2-0.3 mag for geometries 
where the dust is mixed with the gas. In this scenario, the true reddening may be much larger than predicted by E(i? — y). 
We will investigate this possibility in a future study of an unbiased sample of galaxy pairs and N-tuples (Barton, Geller, 
& Kenyon 2000). 

6.2.6. Metallicity 

To investigate the effect of metallicity, we begin with the mean intrinsic [O II] /Ha ratios: 0.34 ± 0.09 and 0.94 ±0.11 
for the H02 and T02 samples, respectively (assiuning an average Ay ^ 1 mag). Using our theoretical grids (equation 12) 
and assuming that the average ionization parameter is g = 3 x 10*" cm/s, the H02 and T02 [O II]/Ha ratios correspond 
to metallicities of log(0/H) -|- 12 ~ 9.18 and 9.06 using our theoretical models, or log(0/H) -|- 12 ~ 8.95 and 8.75 using 
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Fig. 18.— Bottom panels: Comparison between the SFR(Ha) and SFR([0 II]) for the T02 (Tresse et al. 2002) (filled) and the H02 (Hicks et 
al. 2002) (unfilled) samples. SFR([0 II]) and SFR(Ha) are corrected for reddening assuming Av ^ 1. In figure (a), SFR([0 II]) is calculated 
using equation 4 which docs not include an explicit correction for abundance, and assumes log(0/H) + 12~ 8.6. In Figure (b), SFR([0 II]) 
is calculated according to equation 9, with mean M91 metallicities of 8.75 (T02) and 8.95 (H02) derived from the mean [O II]/Hq ratios of 
the two samples. The dotted line shows the least squares fit to the T02 and H02 data. The solid line shows equality between SFR(Ha) and 
SFR([0 II]). Errors in SFR(Ha) and SFR([0 II]) are ~ 30% and ~ 35% respectively, not including errors systematic introduced by using 
Av ~ 1 rather than the true (unknown) extinction. Top panel: The SFR(Ha) versus the logarithm of the ratio of SFR([0 II]) and SFR(Ha) 
from the bottom panel. 

the 1VI91 diagnostic and equation 7. These values are consistent with the recent result by LCS that > 75% of 65 CFRS 

galaxies between 0.5 < z < 1.0 have an oxygen abundance log(0/H) + 12 8.9 using the M91 method. The SFR([0 II]) 
calibration used in Figure 18a assumes that the mean abundance is the same as that measured for the NFGS. The 
absolute value of the abundance is diagnostic dependent (compare the x-axes in Figures 6a - 6d). Because LCS use the 
M91 method, we also use this method to derive the mean NFGS abundance for comparison: log(0/H) + 12 ~ 8.63. The 
difference between the NFGS, T02 and H02 metallicities is likely to result from the different luminosity ranges covered 
by each sample. Most high-z samples consist of more intrinsically luminous galaxies than in local samples. Indeed, when 
LCS compare the mean abundances for the CFRS sample with the NFGS selected over the same luminosity range, they 
find that the mean abundances are similar. 

Hydrodynamic cosmological models yield dilferent predictions about the stellar metallicity distribution as a function 
of redshift. Models by Edvardsson et al. (1993) and Nagamine et al. (2001) predict that the average stellar metallicity 
is constant up to z ^ 2; the models by Rocha-Pinto et al. (2000) predict that earlier stars are more metal-poor. Either 
way, it is likely that most samples observed at high-z contain more intrinsically luminous (and thus higher metallicity) 
galaxies than local samples. 

To check whether the high metallicities result from a luminosity selection effect, we compare the mean Mb and metallicity 
with the metallicity-luminosity relations presented in Kobulnicky et al. (2003). The mean Mb (^ ^21) and metallicity 
(log(0/H) + 12 ~ 8.75) for the T02 sample are consistent with the metallicity-luminosity relation for the DGSS 0.6 < Z < 
0.8 galaxies in Kobulnicky et al. On the other hand, the H02 mean Mb (~ —21 and metallicity (log(0/H) 4- 12 ~ 8.95) 
are more similar to the upper end of the luminosity-metallicity relation in local samples. This result may be caused by 
the H02 Ha sample selection. Such a selection could potentially bias the sample towards galaxies with intrinsically strong 
Ha fluxes, small [O II] /Ha ratios and high metallicities. 

We obtain a metallicity-corrected SFR([0 II]) with either the mean [O II] /Ha ratios and equation (1), or equivalently, 
we use the mean metallicities derived above and equation (15) or (10). Figure (18b) compares the metallicity-corrected 
SFR([0 II]) estimates to SFR(Ha). The SFR([0 II]) and SFR(Ha) agree with a mean relative scatter of 0.17 dex. 
We conclude that correcting for a mean metallicity consistent with the mean [O II] /Ha for the T02 and H02 samples 
produces better agreement between SFR([0 II]) and SFR(Ha). The scatter is identical to that observed for the NFGS 
with Ay ^ 1 mag (Figure 14). 
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7. SFR([0 II]) AND THE COSMIC STAR FORMATION HISTORY 

Analysis of the star formation history of the Universe depends on use of the [O II] emission-line as a SFR diagnostic in 
the z - 0.4 - 1.6 range (e.g., Hammer ct al. 1997; Hogg ot al. 1998; Sullivan ct al. 2000; Gallego et al. 2002; Hippelein 
et al. 2003; Teplitz et al. 2003). There are many approaches to calculating the star formation history based on the [O II] 
luminosity. One approach (Hogg et al. 1998) is to calculate the [O II] luminosity density as a function of redshift and 
then to convert this estimate directly into a star formation rate based on one of the SFR([0 II]) conversions. Hogg et al. 
use the K92 calibration, and conclude that it (the K92 conversion) may not be directly applicable to galaxies at high- 2:. 

Hippcloin ot al. (2003) take another approach. Hippelein ct al. use an extinction-corrected [O II] /Hq: flux ratio (0.9) 
to convert the K98 Ha SFR calibration into an [O II] SFR calibration for their sample which covers the redshift range 
0.25 < z < 1.2. The resulting constant used to convert the intrinsic [O II] luminosity (in ergs/s) into a SFR is 8.8 x 10~^^. 
If the [O II] /Ha ratio that Hippelein et al. derive is a good approximation to the average for their entire sample, then 
the constant 8.8 x 10^^^ may be used to predict a mean abimdance. Equation (10) gives a mean abundance of ~ 8.87 
for the Hippelein et al. sample with the M91 abundance diagnostic. This abimdance estimate is similar to the mean 
abundance of the LCS sample (~ 8.9), and is reasonable if the Hippelein et al. sample contains relatively more luminous 
galaxies than in local samples. An alternative explanation of the low [O II] /Ha ratio observed by Hippelein et al. is that 
the ionization state of the gas dominates the [O II] /Ha ratio and that the average ionization parameter is higher in the 
Hippelein et al. sample than in the NFGS. 

Teplitz et al. (2003) followed a similar method for their sample of 71 0.46 < z < 1.415 galaxies from the STIS Parallel 
Survey. They used the Jansen, Franx, & Fabricant (2001) NFGS relation between absolute blue magnitude and [O II]/Ha 
to derive an [O II] /Ha ratio (uncorrected for reddening) for various luminosity ranges in their sample. Teplitz et al. note 
that the [O II] /Ha ratio is highly dependent on the metallicity and reddening of each individual galaxy. Teplitz et al. 
calculate a mean uncorrected [O II]/Ha ratio of 0.45, lower than those observed locally (Figure 3a). From these ratios 
they convert the [O II] luminosity into an Ha luminosity using the K98 SFR(Ha) relation. Teplitz et al. show that 
the comoving star formation density estimated using Ha systematically exceeds the [O II] estimate derived for the same 
redshift range. However, Teplitz et al. use SFR densities which have been corrected for reddening in various ways. They 
use the uncorrected [O II] luminosity densities with the empirically derived [O II] /Ha to convert the [O II] luminosity 
densities into Ha luminosity densities. This process is similar to the one used by K98 to derive his SFR([0 II]) indicator: 
the difference in reddening at [O II] and at Ha is taken into account, but the reddening of the Ha luminosity is not. 
The derived luminosity densities should therefore be corrected for reddening to allow comparison of the resulting SFR 
densities with those calculated from Ha surveys. Typically, Ha surveys are already corrected for reddening based on some 
assumption about the average attenuation. Many studies of SFR densities based on Ha make a correction for reddening 
based on an assumed Ay ~ 1 mag. 

Figure (19a) shows the discrepancy observed by Teplitz et al. (2003) between star formation densities measured with 
Ha and those measured with [O II]. All data on Figure (19) have been converted to the standard cosmology {h = 0.72, 
= 0.29; Table 6). Some of the SFR densities from Figure (19a) have been corrected or partially corrected for reddening; 
others are uncorrected. 

Gallego et al. (2002) emphasize the importance of using the same assumptions for star formation rate conversion and 
reddening for all data points in SFR density comparisons. For example, Figure (5) in Teplitz et al. (2003) and Figure (19a) 
show that the Ha SFR density point of Pascual et al. (2001) significantly exceeds that of Tresse & Maddox (1998). The 
source of this apparent difference is in the assumed reddening. Pascual et al. (2001) use an Aua ~ 1 which corresponds 
to an E(B-V)~ 0.6 and an attenuation factor of ioO-iKH<,)E{B-v) ^ 3 7. j^^ssq et al assume an average Ay ~ 1 which 
corresponds to an E(B-V)'^ 0.3, corresponding to an average attenuation of ^ 2.0. If we assume the same reddening of 
Ay ~ 1 for both samples, the two Ha SFR density estimates agree to within 10% (see Figure 19b). 

Figure (19a) also shows [O II] SFR densities, including those by Hammer et al. (1997). Because Hammer et al. 
provide the comoving luminosity density (uncorrected for reddening) as a function of redshift, we convert these into 
reddening corrected SFR densities using our SFR([0 II]) formula (equation 4), assuming an average reddening of Ay ~ 1 
. Figure (19b) shows the new SFR([0 II]) densities. The resulting intrinsic SFR([0 II]) densities are ^ 50% larger than 
the SFR([0 II]) densities in Figure (19a), bringing the SFR([0 II]) values into closer agreement with the SFR(Ha) data 
points. This difference is a result of reddening correction. This conclusion can easily be verified by correcting the SFR(Ha) 
densities (calculated by Teplitz et al. using the Hammer et al. data and [O II]/Ha= 0.45) for reddening at Ha using the 
same extinction curve and assuming an Ay ~ 1. Even with a 50% increase, the [O II] SFR density estimates are still 
lower than the Ha SFR estimates, reminiscent of the T02 and H02 galaxies in Section 6. 

So far, we have used equation (4) to estimate the SFR([0 II]) densities. The use of this equation assumes that the 
average abundance for the samples at high- 2; is log(0/H) + 12~ 8.6 (M91 diagnostic) as observed in the NFGS. The LCS 
sample, and most other high-z samples, contain many intrinsically more luminous galaxies than in the NFGS. Galaxies 
with luminosities representative of the local luminosity function are often too faint to be included in high redshift samples. 
As we have disciissed, the star formation rate density for any particular redshift is estimated using the [O II] luminosity 
density and some [O II] SFR calibration. The use of any SFR [O II] calibration requires the calculation of or an assumption 
about the [O II] /Ha ratio. In SFR history studies, the luminosity density is, in principle, corrected for the missing lower 
luminosity galaxies but the assumed [O II] /Ha is not corrected. The [O II] /Ha ratio for high-z samples will be typical 
of the high luminosity (high-metallicity) galaxies observable, despite the fact that the mean metallicity for high redshift 
galaxies must actually be lower than is observed locally. Evidence for an [O II] /Ha ratio typical of high luminosity galaxies 
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Fig. 19. — Top (a): The cosmic star formation history where the SFR densities arc calculated according to Tcplitz et al. (2003); SFR(Ha) 
density is calculated using the K98 SFR(Ha) cahbration, SFR([0 II]) density is calculated using [OII]/Hq = 0.45 and the K98 SFR(Ha) 
calibration. Reddening is the same as used in Tcplitz et al. Middle (b): The SFR densities calculated using a consistent reddening correction. 
Wc used the K98 SFR(Ha) calibration, assuming an A(v)=l (equation 4;Z=log(0/H)+12 = 8.6). Bottom (c): The SFR densities calculated 
using a consistent reddening and metallicity correction; Z=log(0/H)+12 = 8.8 (equation 19). References correspond to: Hammer et al. 

(1997) (H97), Teplitz et al. (2003) (T03), Yan et al. (1999) (Y99), Tresse et al. (2002) (T02), Pascual et al. (2001) (POl), Tresse & Maddox 

(1998) (T98). 
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is easily observed in the Hippclcin ct al. and Tcplitz et al. samples. The mean reddening-corrected [O II] /Ha ratio for 
the Hippelein sample (0.9) corresponds to a mean abundance of log(0/H) + 12^ 8.87 (M91 method). The Teplitz sample 
has a mean uncorrected [O II]/Ha ratio of ~ 0.45. If we assume an Ay ~ 1 (as in Tresse et al. 2002; Pascual et al. 2001), 
then the mean reddening-corrected [O II] /Ha ratio is ~ 0.83. This [O II] /Ha ratio corresponds to a mean metallicity of 
log(0/H) + 12^ 8.90, similar to the Hippelein sample. 

Both Teplitz et al. and Hippelein et al. use an [O II] /Ha ratio typical for high-z galaxies in their SFR density calculation. 
This [O II] /Ha ratio effectively takes the average metallicity into account. We should therefore apply a similar correction. 
If we assume that the average oxygen abundance for the Hammer et al. and Teplitz et al. galaxies is log(0/H) + 12~ 8.9 
(using the M91 method) then equation 4 becomes: 

SFR(M0yr-i) = (9.53 ± 0.91) x lO"''^ l([OII]) (ergs s'^) (19) 

where we have assumed the NFGS rms dispersion about the mean abundance. The L([0 II]) is corrected for reddening 
assuming Ay ~ 1 and a foreground screen geometry. The resulting SFR densities arc a factor of 1.45 larger than the 
estimates with equation (4). The [O II] SFR density estimates are still slightly lower than the Ha SFR density estimates, 
but the the Ha and [O II] SFR densities now agree to within ~ 30% (Figure 19c). 

An alternative explanation for the difference between the [O II] and Ha SFR densities in Figure 19b is that the dust 
geometry may not be a simple screen. Witt & Gordon (2000) show that E{B — V) saturates at around 0.2-0.3 mag for 
geometries where the dust is mixed with the gas. In this scenario, the true reddening may be much larger than predicted 
by E{B - V). 

We conclude that the major differences between Ha and [O II]-based estimates of the star formation rate density as a 
function of redshift probably result from inconsistent assumptions about reddening and abundance, and/or from failure to 
correct for the biases they introduce. Careful correction for both the reddening and abundance brings the star formation 
rate density estimates into agreement to within ^ 30%. 

8. CONCLUSIONS 

We investigate the use of the [O II] emission-line as a star formation indicator for a sample of 97 nearby field galaxies. 
Our high S/N integrated (global) spectra allow correction of the Balmer lines for underlying stellar absorption and 
correction of the emission-line fluxes for abundance. We find: 

• There is a systematic difference between the SFR(Ha) and SFR([0 II]) using the Kennicutt (1998) calibrations. 
The difference results from the use of the observed (uncorrected) [O II] /Ha ratio to convert the [O II] luminosity 
into a SFR indicator. 

• We derive a new SFR([0 II]) indicator (equation 4) which is independent of the reddening between [O II] and Ha. 
This indicator removes the systematic difference observed for the Kennicutt calibration. 

We estimate the abundances for our sample using four independent abundance diagnostics: Kewley & Dopita (2002) 
([N II]/[0 II]), McGaugh (1991) (R23), Zaritsky, Kennicutt & Huchra (1994) (R23), and Chariot & Longhetti (2001) 
([O IIl]/H/3). We show that: 

1. There is a strong correlation between abundance and the [O II] /Ha ratio. 

2. There is no strong dependence of [O II] /Ha on the ionization parameter for log(0/H) + 12> 8.5 

3. The observed discrepancy between the Kennicutt (1998) SFR(Ha) and SFR([0 II]) calibrations result from redden- 
ing and abundance differences between the K92 sample and other surveys. 

Our theoretical stellar population synthesis and photoionization models support these conclusions. We derive a new 
SFR([0 II], Z) calibration which includes a correction for abundance (equation 10). This relation reduces the rms residuals 
between SFR(Ha) and SFR([0 II]) significantly from 0.11 to 0.03-0.05 dex. Our SFR([0 Il],Z) calibration is based on the 
intrinsic (reddening corrected) [O II] luminosity and the abundance. This calibration can be used to derive a SFR relation 
for samples where the mean abundance differs from the NFGS. The SFR([0 II], Z) relation may also be applied to samples 
where abundance estimates are available, but where H/3 has too low a S/N to be used to calculate a star formation rate 
based on Ha. 

We observe a strong correlation between the intrinsic [O II] luminosity and E(B-V) for the NFGS. Using the observed 
[O II] luminosity to estimate the reddening produces no systematic shift between the SFR(Ha) and SFR([0 II]), but there 
is a large dispersion (0.21 dex). Further investigation is required to verify whether a similar reddening-luminosity relation 
holds for other samples. 

We discuss the application of our SFR([0 II]) calibrations to studies of the cosmic star formation history. We show 
that we can remove the previously observed discrepancy between the star formation rates based on [O II] and those based 
on Ha by correcting for reddening and abundance. The agreement between SFR(Ha) and SFR([0 II], Z) emphasizes the 
importance of aquiring data which yield internally consistent reddening and abundance estimates over the entire redshift 
range of the sample. 
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27 


31 
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181 


NGC 7077 


1142 


0.15 


0.13 


8.55 


8,56 


8,52 


8,56 


6.50 


6.62 


0.10 


0.16 


0.10 


0.08 


0.07 


183 


A22306+0750 


1995 


0.57 


0.16 


8.87 


8.95 


8.77 


8.99 


7.46 


7.51 


0.86 


0.54 


0.82 


0.99 


0.20 


186 


A22426+0610 


1925 


0.16 


0.25 


8.81 


8.92 


8.78 


8.71 


7.17 


7.08 


0.44 


0.46 


0.30 


0.34 


0.45 


187 


A22551 + 1931N 


5682 


0.61 


0.30 


8.90 


8.93 


8.77 


8.80 


8.12 


8.12 


4.01 


2.00 


3.32 


3.79 


1.12 


189 


NGC 7460 


3296 


0.48 


0.31 


8.97 


9,10 


8.91 


9.13 


8.01 


7.90 


3.12 


1.58 


1.98 


4.14 


1.56 


190 


NGC 7537 


2648 


0.41 


0.31 


8.93 


8,96 


8,79 


8.87 


7.79 


7.80 


1.87 


1.43 


1.57 


1.94 


1.23 


192 


A23176+1541 


4380 


0.49 


0.34 


8.81 


8.84 


8.69 


8.81 


7.95 


8.06 


2.71 


2.24 


2.87 


2.74 


1.56 


193 


NGC 7620 


9565 


0.56 


0.59 


8.88 


8.92 


8.76 


8.85 


9.11 


9.14 


38.98 


23.30 


34.87 


39.65 


41.12 


195 


IC 1504 


6306 


0.92 


0.35 


8.78 


8.77 


8.57 


9.04 


8.69 


8.92 


14.86 


6.76 


20.94 


18.08 


2.05 


196 


NGC 7752 


4902 


0.37 


0.42 


8.97 


9.04 


8.87 


8.84 


8.25 


8.12 


5.43 


3.31 


3.32 


5.30 


5.71 


198 


A23542+1633 


1788 


0.11 


0.26 


8.11 




8.08 


8.56 


6.97 


7.03 


0.28 


0.46 


0.27 


0.20 


0.38 



^E(B-V) is the observed reddening calculated using the Balmer Decrement as described in Section 2. E(B-V) includes Galactic extinction which is ~ 0.015 on average for the NFGS 
^E(B-V) is calculated using the uncorrected (observed) [Oil] luminosity from equations (16) and (18), described in Section 6. 

'^Thc Hq luminosity is corrected for underlying stellar absorption and for reddening using the E(B-V) derived from the Balmer decrement. Units are in the logarithmic scale in term of Lq 
'^The [on] luminosity is corrected for reddening using the E(B-V) derived from the Balmer decrement. Units are in the logarithmic scale in terms of Lq 
^SFR([OII]) is derived using equation (4) and L([OII]) corrected for reddening using the E(B-V) derived from the Balmer decrement. 

^SFR([OII]) is derived using equation (15) with the Z94 abundance. R23 and L([OII]) have been corrected for reddening using the E(B-V) derived from the Balmer decrement. 
®SFR([OII]) is derived using the intrinsic luminosity estimated from the observed luminosity (equation 18). Equation (10) was used to derive the SFR([OII]). 
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Table 2 

Slope, y-intercept, and scatter for the line of best fit in equation 7 (Figures 6a-d). 



Figure 


Abundance 


slope 


y-intercept 


rms 


Spearman Rank 


Number 


Diagnostic 


a 


b 


a 


Coefficient Probability (%) 


6a 


KD02 [N II]/[0 II] 


-1.99 ±0.32 


18.67 ±2.83 


0.10 


-0.79 1.75 X lO"^*^ 


6b 


Z94 R23 


-1.75 ±0.25 


16.73 ±2.23 


0.07 


-0.88 3.80 X lO"^** 


6c 


M91 R23 


-2.29 ±0.41 


21.21 ± 3.66 


0.04 


-0.93 < 10"^° 


6d 


COl [O III]/H/3 


-0.88 ±0.17 


8.78 ± 1.53 


0.24 


-0.36 0.05 



Table 3 

Slope, y-intercept, and scatter for the line of best fit to the SFR(Ha) versus SFR([0 II]) plots. 



Figure 


Abundance 


slope 


y-intercept 


rms 


Number 


Diagnostic 





b 


a 


2 (K98) 


none 


0.83 ±0.02 


0.01 ±0.02 


0.11 


4 


none 


0.97 ±0.02 


-0.03 ±0.02 


0.08 


7a 


KD02 IN II] /[O II] 


1.04 ±0.03 


-0.007 ±0.03 


0.05 


7b 


Z94 R23 


1.02 ±0.03 


0.01 ±0.03 


0.04 


7c 


M91 R23 


1.01 ±0.03 


-0.012 ±0.03 


0.03 


7d 


COl [O III]/H/3 


1.01 ±0.03 


0.04 ± 0.03 


0.08 


12a 


KD02 [N II] /[O II] 


1.03 ± 0.04 


-0.04 ± 0.03 


0.05 


12b 


Z94 R23 


1.03 ±0.02 


+0.03 ± 0.01 


0.05 


14a 


Z94 (Av = 1) 


0.76 ± 0.03 


-0.06 ±0.02 


0.17 


14b 


Grid (Av = 1) 


0.77 ±0.03 


-0.04 ±0.02 


0.17 


16a,b 


Z94; E(B-V) from L(JO II]) 


1.07 ±0.03 


-0.06 ±0.03 


0.21 


18a 


none (Av = 1) 


0.93 ±0.12 


-0.20 ±0.11 


0.21 


18b 


M91 mean (Ay = 1) 


1.09 ±0.11 


-0.15 ±0.10 


0.17 
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Table 4 

Coefficients for the best fit to the theoretical curves for abundance versus [O II] /Ha. 



q (cm/s) 


a 


6 


c 


d 


5x10® 


-2564.67 


847.554 


-92.9404 


3.38261 


1x10'^ 


-2877.94 


955.234 


-105.245 


3.85016 


2x10^ 


-2281.80 


754.840 


-82.8388 


3.01667 


4x10'^ 


-1432.67 


470.545 


-51.2139 


1.84750 


8x10'^ 


-786.096 


256.059 


-27.6048 


0.984957 



Table 5 



Luminosity and SFRs for the seven Hicks et al. (2002; H02) galaxies and 30 Tresse et al. (2002; T02) galaxies 

WITH [O II] AND Ha luminosities ( Av = 1 assumed) . 



ID 


Ref 


z 


log 


F[OII] 

(10~^*ergs" 


F(Ha) 
■1 cm-2) 


log 


log 


log 


log 


SFR[OII] 


SFR[OII] 


SFR 








L([OII])o 


L([OII])* 


L(Ha)o 


L(Ha)i 


(8.6)» 


(8.75,8.95)'= 


Ha 


2 


H02 


0.76 


-0.56 


98.61 


358.03 


41 


.40 


41 


.97 


41.96 


42.27 


6.20 


21.91 


14.58 


4 


H02 


1.14 


-0.39 


79.21 


194.44 


41 


.74 


42 


.32 


42.13 


42.44 


13.60 


48.02 


21.61 


5 


H02 


1.45 


-1.07 


5.64 


66.26 


40 


.85 


41 


.43 


41.92 


42.23 


1.76 


6.21 


13.38 


7 


H02 


1.06 


-0.71 


23.27 


119.34 


41 


.13 


41 


.70 


41.84 


42.15 


3.33 


11.78 


11.07 


9 


H02 


1.49 


-1.05 


10.67 


119.72 


41 


.16 


41 


.73 


42.21 


42.51 


3.56 


12.57 


25.85 


10 


H02 


1.37 


-1.03 


11.81 


126.55 


41, 


.11 


41, 


.69 


42.14 


42.45 


3.20 


11.30 


22.19 


11 


H02 


1.13 


-0.90 


11.81 


93.81 


40 


.91 


41 


.48 


41.81 


42.11 


1.98 


7.01 


10.20 


0.0338 


T02 


1.03 


-0.20 


113.90 


178.50 


41 


.79 


42 


.36 


41.98 


42.29 


15.12 


19.32 


15.34 


0.0564 


T02 


0.61 


-0.67 


30.10 


140.30 


40 


.65 


41 


.23 


41.32 


41.63 


1.11 


1.41 


3.34 


0.0874 


T02 


0.71 


-0.46 


22.30 


64.50 


40 


.69 


41 


.26 


41.15 


41.45 


1.19 


1.53 


2.24 


3.0167 


T02 


0.60 


-0.09 


57.30 


71.30 


40, 


.92 


41, 


.49 


41.01 


41.32 


2.04 


2.60 


1.64 


3.0422 


T02 


0.71 


-0.72 


50.10 


265.60 


41 


.04 


41 


.61 


41.77 


42.07 


2.71 


3.46 


9.30 


3.0480 


T02 


0.61 


-0.19 


155.60 


238.50 


41 


.36 


41 


.93 


41.55 


41.85 


5.66 


7.24 


5.62 


3.0485 


T02 


0.61 


-0.52 


249.70 


823.20 


41 


.56 


42 


.14 


42.08 


42.39 


9.01 


11.50 


19.22 


3.0488 


T02 


0.61 


-0.21 


278.20 


450.50 


41 


.61 


42 


.19 


41.82 


42.13 


10.09 


12.88 


10.57 


3.0570 


T02 


0.65 


-0.37 


94.00 


220.70 


41, 


.21 


41, 


.78 


41.58 


41.89 


4.00 


5.11 


6.08 


3.0589 


T02 


0.72 


-0.60 


57.00 


228.90 


41 


.10 


41 


.67 


41.70 


42.01 


3.10 


3.96 


8.05 


3.0595 


T02 


0.61 


-0.69 


54.40 


269.40 


40 


.90 


41 


.48 


41.60 


41.90 


1.97 


2.51 


6.30 


3.0615 


T02 


1.05 


-0.20 


88.10 


139.10 


41 


.70 


42 


.27 


41.90 


42.20 


12.27 


15.68 


12.54 


3.0879 


T02 


0.60 


-0.16 


41.10 


59.40 


40 


.77 


41 


.35 


40.93 


41.24 


1.46 


1.86 


1.36 


3.0952 


T02 


0.86 


-0.12 


58.10 


76.00 


41 


.30 


41 


.87 


41.42 


41.72 


4.93 


6.30 


4.18 


3.0984 


T02 


0.70 


-0.85 


20.40 


144.50 


40 


.63 


41 


.21 


41.48 


41.79 


1.06 


1.35 


4.84 


3.1027 


T02 


1.04 


-0.19 


106.20 


164.40 


41 


.77 


42 


.34 


41.96 


42.26 


14.45 


18.45 


14.48 


3.1032 


T02 


0.62 


0.28 


84.10 


44.60 


41 


.11 


41 


.69 


40.84 


41.14 


3.19 


4.07 


1.09 


3.1242 


T02 


0.77 


-0.51 


66.90 


214.30 


41 


.25 


41 


.82 


41.75 


42.06 


4.33 


5.53 


8.98 


3.1309 


T02 


0.62 


-0.54 


104.10 


364.10 


41 


.20 


41 


.78 


41.75 


42.05 


3.93 


5.02 


8.90 


3.1345 


T02 


0.62 


-0.24 


70.00 


121.90 


41 


.03 


41 


.60 


41.27 


41.58 


2.64 


3.37 


2.97 


3.1349 


T02 


0.62 


-0.30 


128.60 


255.20 


41 


.29 


41 


.87 


41.59 


41.89 


4.83 


6.16 


6.20 


3.1534 


T02 


0.80 


-0.00 


127.90 


128.10 


41 


.57 


42 


.14 


41.57 


41.87 


9.07 


11.58 


5.88 


3.9003 


T02 


0.62 


-0.48 


187.70 


564.10 


41 


.46 


42 


.04 


41.94 


42.24 


7.14 


9.12 


13.89 


22.0429 


T02 


0.63 


-0.60 


65.30 


261.20 


41 


.02 


41 


.59 


41.62 


41.92 


2.56 


3.27 


6.63 


22.0637 


T02 


0.54 


-0.56 


151.60 


554.80 


41 


.24 


41 


.81 


41.80 


42.10 


4.23 


5.40 


10.02 


22.0764 


T02 


0.82 


-0.36 


51.30 


118.70 


41, 


.20 


41 


.77 


41.56 


41.87 


3.89 


4.96 


5.82 


22.0770 


T02 


0.82 


-0.31 


100.40 


206.70 


41 


.49 


42 


.06 


41.80 


42.11 


7.59 


9.70 


10.12 


22.0779 


T02 


0.93 


-0.26 


38.60 


69.60 


41 


.21 


41 


.78 


41.46 


41.77 


3.95 


5.04 


4.61 


22.0843 


T02 


0.92 


-0.66 


37.00 


170.60 


41 


.18 


41 


.75 


41.84 


42.14 


3.68 


4.70 


10.99 


22.1406 


T02 


0.82 


-0.28 


175.10 


333.50 


41 


.73 


42 


.30 


42.01 


42.31 


13.21 


16.88 


16.29 



^SFR([OII]) is calculated using equation (4). This equation assumes a metallicity of log(0/H) + 12 = 8.6 (the mean metallicity of our NFGS sample 
using the M91 method). 

^SFR([OII]) is calculated using equation (9) with log(0/H) + 12 = 8.75 (M91 method) for the T02 sample and log(0/H) + 12 = 8.95 (M91 method) 
for the H02 sample 



34 



Table 6 

StAH FOKM.VTIOX mate DEKSITIES IK FKU'WE 19. 



z 


Dataset 


line 


SFR density (Mq yr'^ Mpc"^) 








orig 


new 
(Z=8.6) 


new'^ 
(Z=8.9) 


0.20 


Tresse & Maddox 1998 (M98) 


Ha 


0.01 


0.01 


0.01 


0.24 


Pascual et al. 2001 (POl) 


Ha 


0.03 


0.01 


0.01 


0.38 


Hammer et al. 1997 (H97) 


[Oil] 


0.005 


0.007 


0.01 


0.60 


Tresse ct al. 2002 (T02) 


Ha 


0.07 


0.07 


0.07 


0.63 


Hammer et al. 1997 (H97) 


[Oil] 


0.02 


0.03 


0.04 


0.72 


Teplitz et al. 2003 (T03) 


[OH] 


0.03 


0.04 


0.06 


0.80 


Tresse et al. 2002 (T02) 


Ha 


0.10 


0.10 


0.10 


0.88 


Hammer et al. 1997 (H97) 


[OH] 


0.05 


0.07 


0.11 


0.94 


Teplitz et al. 2003 (T03) 


[OH] 


0.04 


0.06 


0.09 


1.14 


Teplitz et al. 2003 (T03) 


[OH] 


0.07 


0.09 


0.14 


1.30 


Yan et al. 1999 (Y99) 


Ha 


0.12 


0.24 


0.24 



'^SFR densities arc calculated according to Teplitz et al. (2003): SFR(Ha) density is calculated using the K98 SFR(Ha) calibration, 
SFR([OII]) density is calculated using [OH] /Ha = 0.45 and the K98 SFR(Ha) calibration. Reddening is kept as in Teplitz et al. 

and is not the same for each data set. 

''SFR(Ha) densities calculated using the K98 SFR(Ha) calibration, assuming an A(v)=l. SFR([OII]) densities calculated using 
our SFR([OII]) calibration (equation 4), assuming A(v)=l and Z=log(0/H)+12 = 8.6. 

'^SFR(Ha) densities calculated using the K98 SFR(Ha) calibration, assuming an A(v)=l. SFR([OII]) densities calculated using 
our SFR([OII]) calibration (equation 19), assuming A(v)=l and Z=log(0/H)+12 = 8.9. 



